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INTRODUCTION 

The  objective  of  the  project  is  to  test  whether  the  therapeutic  effectiveness  of 
permanent  implant  brachytherapy  for  prostate  cancer  can  be  improved  by  using  a 
combination  of  short  and  long  half  life  radionuclides  simultaneously.  Specific  aims  of 
the  proposed  project  are: 

1.  To  test  theoretically  the  potential  of  a  mixture  of  radionuclides  in  permanent 
implants,  using  the  linear  quadratic  model,  as  a  function  of  Tpot,  potential  tumor  doubling 
time. 

2.  To  test  experimentally  the  validity  of  this  concept  by  in  vitro  irradiation  of 
BA1 112  sarcoma  cells  at  a  continuous  low  dose  rate  (CLDR)  with  125I  (60  d  half  life), 
103Pd  (17  d  half  life)  and  a  50:50  mixture  of  12'T  and  103Pd  under  aerobic  conditions 
leading  to  exponential  growth  at  different  rates  (from  near  quiescence  to  full  exponential 
growth  at  a  maximal  rate,  with  a  doubling  time  of  approximately  14  hours). 

3.  To  measure  the  radiobiology  parameters  such  as  alpha,  beta,  half  life  of  repair  for 
the  BA1 1 12  sarcoma  cells  under  different  growth  conditions  and  develop  a  theoretical 
model  to  predict  expected  levels  of  cell  killing  using  1251,  103Pd  or  a  mixture  of  these 
isotopes. 

4.  To  use  immunohistochemical  techniques  to  measure,  in  solid  BA1 1 12  tumors  in 
vivo,  the  proportion  of  cells  in  S  phase,  the  proportion  proliferating  and  non-proliferating 
cells  and  the  tumor  doubling  time. 

5.  To  test  the  therapeutic  effectiveness  of  103Pd,  125I  and  a  Pd/I  mixture  in  the 
BA1 1 12  in  vivo  tumor  system; 

6.  To  test  the  therapeutic  effectiveness  of  103Pd,  125I  and  a  Pd/I  mixture  in  human 
prostate  carcinoma  xenografts  in  nude  mice,  using  a  slow  growing  and  a  fast  growing 
carcinoma. 

7.  To  evaluate  the  clinical  potential  and  feasibility  of  this  approach  in  the  treatment 
of  human  prostate  cancer. 

BODY  OF  THE  REPORT 

We  have  developed  a  theoretical  radiobiology  model  for  cell  killing  by 
continuous  low  dose  rate  irradiation  (CLDRI)  using  a  mixture  of  radionuclides. 
Theoretical  studies  were  performed  to  investigate  the  hypothesis  and  to  plan  in  vitro  and 
animal  studies.  Experiments  have  been  performed  using  BA1112  tumor  cells  and  Chinese 
Hamster  cells  growing  in  vitro  and  BA1112  cells  growing  in  vivo  as  solid  tumors  in 
WAG/rij  rats.  Radiobiology  parameters  for  these  cells  have  been  determined  and  used  in 
the  theoretical  radiobiology  model  to  improve  our  understanding  of  the  experimental 
observations. 
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The  linear-quadratic  model  of  cell-killing  by  CLDBI  Using  a  Mixture  of 
Radionuclides 


We  have  developed  a  theoretical  model  for  CLDRI  using  a  mixture  of 
radionuclides  with  different  half  lives.  This  model  is  described  in  the  attached  manuscript 
entitled  “Biologically  effective  dose  (BED)  for  interstitial  seed  implants  containing  a 
mixture  of  radio-nuclides  with  different  half  lives”  by  Zhe  Chen  and  Ravinder  Nath. 
Briefly,  the  purpose  of  this  project  was  to  develop  a  tool  for  evaluating  interstitial  seed 
implants  that  contain  a  mixture  of  radionuclides  with  different  half-lives  and  to  examine 
the  clinical  implications  of  prescribing  to  an  isodose  surface  for  such  an  implant.  Using  a 
generalized  equation  for  the  biological  effective  dose  (BED)14,  the  effects  of  cell 
proliferation  and  sub-lethal  damage  repair  were  examined  systematically  for  implants 

containing  a  mixture 
of  radionuclides 
(Figure  1).  The 
results  were 
contrasted  with  those 
for  implants  using  a 
single  type  of 
radionuclide.  A 
clinical  permanent 
seed  implant  that 
contained  a  mixture 
of  125I  and  103Pd  seeds 
was  used  to  examine 
the  clinical 
implications  of  the 
isodose  prescription 
for  such  implants.  An 
equation  of  BED  for 
implants  containing 
any  number  of 
radionuclide  types 

was  obtained.  For  implants  containing  a  mixture  of  radionuclides  with  different  half-lives 
such  as  125I  and  103Pd,  the  dose  as  well  as  its  temporal  delivery  pattern  to  a  point  is 
dependent  on  the  relative  dose  contributions  from  different  types  of  radionuclide.  It  can 
vary  from  point  to  point  throughout  the  implant  volume.  Therefore  the  quantitative 
effects  of  cell  proliferation  and  sub-lethal  damage  repair  are  spatially  dependent  in  such 
an  implant.  For  implants  containing  a  mixture  of  125I  and  103Pd  seeds,  the  prescription  to 
an  isodose  surface  becomes  non-unique.  If  the  prescription  dose  was  based  on  existing 
clinical  experience  of  using  125I  seeds  alone,  mixing  lfePd  seeds  with  125I  seeds  would 
decrease  the  cell  survival  in  such  implant.  On  the  other  hand,  if  the  prescription  dose  was 
based  on  existing  clinical  experience  of  using  103Pd  seeds  alone,  mixing  12U  seeds  with 
l03Pd  seeds  in  a  same  implant  would  create  radiobiologically  “cold”  spots  (i.e.  an  increase 
in  cell  survival  from  the  clinical  expectation)  at  locations  where  a  major  portion  of 
prescription  dose  is  contributed  by  the  125I  seeds.  For  fast-growing  tumors,  these  "cold" 
spots  can  become  significant.  From  this  theoretical  investigation,  we  conclude  that  when 
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cell  proliferation  and  sub-lethal  damage  repair  are  present  during  dose  delivery,  total  dose 
alone  is  no  longer  sufficient  for  a  complete  characterization  of  an  interstitial  seed  implant. 
In  order  to  avoid  radiobiological  “cold”  spots  when  radionuclides  of  different  half-lives 
are  mixed  in  a  permanent  implant,  the  dose  prescription  should  be  based  on  the  clinical 
experience  of  using  the  longer  half-life  radionuclide.  Biologically  effective  dose  provides 
a  tool  to  start  examining  the  radiobiological  effects  of  mixing  different  type  of 
radionuclides  in  the  same  implant. 

A  manuscript  has  been  submitted  for  publication  in  the  International  Journal  of  Radiation 
Oncology. 


In  vitro  CLDRI  studies 

BA1 1 12  tumors  were  grown  between  the  ears  of  the  a  14  week  old  male  WAG/rij 
rats  by  interdermal  inoculation  from  a  single  cell  suspension  of  B All  12  cells  obtained 
from  a  21  day  BA1 112  tumor  growing  on  the  head  of  a  previously  inoculated  rat.  A 
tumor  cell  suspension  was  made  from  the  BA1 1 12  tumor  and  between  1.5x10s  and 
5.0x10s  cells  were  plated  into  petri  dishes.  These  cells  were  allowed  to  settle  and  reach 
logarithmic  growth,  (48-72  hrs),  before  they  were  used  in  a  continuous  low  dose  rate 
experiment. 


Monolayers  of  rat 
rhabdomyosarcoma  cells  (BA1112) 
were  irradiated  in  vitro  by  103Pd 
sources  in  a  polystyrene  phantom. 
Colony  formation  ability  of  irradiated 
cells  under  aerobic  conditions  was 
measured  for  graded  doses,  at  a  dose 
rate  of  6  cGy/hr.  Dose  to  the  cell 
monolayers  was  determined  using 
FeSC>4  Fricke  dosimetry,  with  a 
calculated  correction  for  interface 
effects  due  to  photoelectric  effect  in 
the  tissue  culture  dishes.  The  sources 
(up  to  60  in  one  experiment)  were 
arranged  in  concentric  circles  in  such 
a  way  as  to  provide  a  dose  uniformity 
of  better  than  ±  5%  across  the  dishes. 
Some  of  the  results  are  shown  in 
Figure  2.  Comparison  of  the  surviving 
fraction  as  a  function  of  dose 
calculated  for  the  BA1 112  cells  to 
that  measured  by  using  CLDRI  103Pd  irradiation  is  also  shown  in  Figure  2.  The  lines 
through  the  data  points  represent  the  calculated  survival  curves  the  symbols  represent  the 


Using  \ 
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measured  data  The  parameters  of  a,  (3,  repair  half-time,  and  tumor  doubling  time  were 
determined  directly  from  the  measurements  performed  on  the  BA1 1 12  cells,  as  described 
later  in  the  report. 

A  manuscript  entitled  “Dose  rate  dependence  of  the  relative  biological  effectiveness  of 
103Pd  for  Continuous  Low  Dose  Rate  Irradiation  of  B A1 1 12  Rhabdomyosarcoma  Cells  in 
vitro  relative  to  Acute  Exposures”  has  been  prepared  for  submission  to  Int.  J.  Radiat. 
Oncology.  (Appendix  II). 


In  vitro  studies  at  an  acute  dose  rate  using  simulated  x-rav  heams 

To  study  the  radiobiological  characteristics  of  the  cells  under  acute  exposure 
condition,  simulated  x-ray  beams  with  average  energies  equivalent  to  that  emitted  by  125I 
(27.2  -  35.49  keV  with  an  average  of  27.4)  and  103Pd  (20  -  22.7  keV  with  an  average  of 
20.5  keV)  were  established  on  a  new  orthovoltage  unit.  The  simulated  beams  not  only 
have  the  average  energies  similar  to  that  given  by  the  radioactive  isotopes  but  also  have  a 
narrow  photon  energy  spectrum.  The  narrow  photon  energy  spectrum  was  achieved  by 
optimizing  the  tube  voltage  (which  determines  the  upper  limit  of  the  produced  photon 
energy)  and  the  added  filtration  (which  filters  out  the  low-energy  Bremsstrahlung 
photons),  following  the  work  of  Muench  et  al5. 

Aluminum  filters  from  Pantak  was  used  to  construct  a  customized  filter  for  the 
DXT  300  unit  that  has  a  desired  filtration  thickness.  A  set  of  aluminum  filters  (with 
thickness  of  0.1  to  1.0  mm)  from  Nuclear  Associates  (AL  Filter  Set  07-430)  were  used  to 
determining  the  half-value-layer  (HVL)  of  a  simulated  beam  using  a  customized 
filtration.  The  thicknesses  of  the  aluminum  sheets  were  measured  by  using  a  Mitutogo 
micrometer  (Serial  #  2032360)  with  accuracy  of  0.001  mm).  The  aluminum  HVL  for  a 
given  beam  was  determined  by  in-air  ionization  chamber  measurement  under  the  narrow 
beam  geometry.  An  air-equivalent  Spokas  chamber  (Exradin,  Model  No:  Al  (0.5  ml,  AE 
plastic))  was  used  to  measure  the  ionization  at  a  fixed  source  to  chamber  distance  (SCD) 
of  50  cm  in  air.  The  ionization  charge  was  measured  by  a  Keithley  electrometer  (model 
356 14E  SN  43075)  with  -300  V  bias  potential.  Due  the  energy  dependence  of  the 
chamber  at  the  low  energies,  the  measured  ionization  were  converted  to  corresponding 
exposure  and  the  HVL  is  then  determined  from  the  relative  exposure  as  a  function  of 
aluminum  filter  thickness.  A  narrow  circular  beam,  with  a  diameter  of  6  cm  at  SCD  of  50 
cm,  was  generated  by  using  a  homemade  lead  collimator  mounted  to  DXT  300’ s 
accessory  mount.  The  aluminum  sheets  added  to  the  beam  for  the  HVL  measurement 
were  taped  to  the  bottom  of  the  lead  collimator. 

The  HVL  as  a  function  of  mono-energetic  photon  beam  energy  for  aluminum  is 
taken  from  Johns  and  Cunningham6.  The  expected  HVL  for  a  simulated  125I  (103Pd)  beam 
with  average  energy  of  27,4  keV  (20.5  keV)  is  1.84  mm  (0.82  mm)  aluminum.  With  the 
expected  HVL  in  mind,  the  tube  kV,  mA  and  the  thickness  of  the  added  filtration  were 
optimized  for  a  simulated  125I  and  a  simulated  103Pd  x-ray  beam.  The  optimum  setting 
determined  for  the  DXT  300  unit  is  summarized  in  the  following  table  I. 
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Table  I.  Radiation  characteristics  of  simulated  x-ray  beams 


Beam 

kV 

niA 

Added  Filter 
(mm  AL) 

Beam  HVL 
(mm  AL) 

Energy 

Homogeneity 

(%) 

Equivalent 

Energy 

(keV) 

<E>  from 
isotope 
(keV) 

125I 

Equivalent 

43 

20 

3.545 

1.851 

86.9 

27.45 

27.4 

'“Pd 

Equivalent 

29 

25 

1.826 

0.82 

88.6 

20.5 

20.5 

In  vitro  studies  for  quiescent  BA1112  cells 


Animals  were  implanted  with 
transplanted  BA1 112 
rhabdomyosarcomas  by 
inoculation,  into  a  subcutaneous 
site  on  the  heads.  The  tumors 
were  allowed  to  grow  for  3 
weeks,  to  an  experimental 
volume  of  approximately  199- 
200mm3.  The  animals  chosen 
for  the  quiescent  cell 
experiments  were  euthanized  by 
anesthetic  overdose  and  the 
tumor  cells  will  then  be 
removed  using  aseptic 
techniques,  A  single-cell 
suspension  of  tumor  cells  was 
suspended,  counted,  and 
assayed  for  viability  using  the 
same  colony  formation  assay 
used  for  cells  in  cultures.  1.5 
xlO5  cells  were  plated  into  a 
flask  with  13  ml  of  DMEM  for 
cell  growth.  These  cells  will 
then  be  passed  twice  a  week  for  approximately  4-8  passages.  The  cells  are  transplanted 
for  2-4  weeks  to  assure  a  homogeneity  population  of  BA1 112  cells. 

When  the  first  line  has  reached  passage  5  the  cells  were  counted  and  1  .OxlO7  plated  into  5 
petri  dishes  for  a  cell  survival  experiments  with  either  normal  or  quiescent  cells.  We 
know  from  our  previous  quiescent  cell  induction  experiments  that  between  Days  5-7  were 
optimal  to  assure  cells  in  a  healthy  quiescent  state.  If  we  waited  past  Day  7  we  will  be  in 
past  are  optimal  window  and  exhibiting  the  cell  death  phase. 


250  kV  X-Ray  of  Normal  &  Quiescent  Cells 


Dose  (Gy) 
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The  quiescent  cells  are  irradiated  on  Day  6  or  7,  and  the  response  of  tumors  to  irradiation 
will  be  assessed  by  the  ability  of  the  BA  1112  tumor  cells  to  form  colonies  in  cell  culture 
after  in  vitro  irradiations.  A  single-cell  suspension  of  the  tumor  cells  will  be  counted,  and 
assayed  for  viability  and  appropriate  cell  numbers  plated  in  each  petri  dish  as  determined 
by  the  specific  dose  given  to  the  BA1 1 12  cells  in  the  petri  dish.  Therefore  multiple 
platings  from  a  single  cell  suspension  are  mandatory  to  ensure  an  adequate  countable 
number  of  colonies.  After  these  dishes  are  allowed  to  grow  for  14  days  in  a  controlled 
environment,  37°C  incubator  with  95%  air/5%  CO2,  they  are  stained  with  crystal  violet. 
Analysis  of  cell  yield  is  performed. 


Radiological  Parameters  for  CCL-16  and  BA1112  cells 

In  order  to  measure  the  radiological  parameters  of  the  cells  used  in  these  studies,  the 
acute  exposure  survival  curves  were  measured  using  250  kV  x-rays.  Split  dose 
experiments  were  conducted  to  measure  the  half  time  of  sublethal  damage  repair.  Results 
are  described  in  Appendix  III  and  IV.  For  BA1 1 12  cells  the  results  are  summarized  in  the 
table  below. 


Radiobiological  Parameters  Determined  from  Acute  Exposure  Experiments 


Parameters  for  Ba-1112 

Beam 

a 

(oy1) 

P 

(Gy2) 

a/p 

(Gy) 

Tumor  Doubling 
Time 
(day) 

Sub-lethal  Damage 
Repair  Half-time  (hr) 

250  kV  old 

0.25 

0.041 

6.1 

3.0 

0.33 

125I  equivalent 

0.26 

0.043 

6.1 

3.0 

0.33 

103Pd  equivalent 

0.29 

0.048 

6.1 

3.0 

0.33 

250  kV  old 

0.23 

0.044 

5.1 

3.0 

0.33 

mI  New 

0.25 

0.044 

5.6 

3.0 

0.33 

103Pd  New 

0.32 

0.044 

7.3 

3.0 

0.33 

IUdR  and  BrdU  Labeling  For  Flow  Cytometric  Analysis 

We  have  developed  protocols  for  in  vitro  IudR  labeling  and  in  vivo  BrdU  labeling  for 
flow cytometric  analysis  of  CCL-1 16  and  BA1 1 12  cells  based  upon  a  protocol  from  Dr. 
Hong  Sun  at  Yale  (Appendix  V).  Several  studies  on  CCL-1 16  and  BA1 1 12  cells  have 
been  performed  and  results  are  being  analyzed.  An  example  is  shown  below. 
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In  vivo  tumor  growth  studies 

In  order  to  produce  a  consistent  dose  distribution  to  irradiate  the  tumors 
transplanted  to  different  animals  and  to  minimize  radiation  exposure  to  personnel 
handling  the  radioactive  seeds,  an  afterloading  seed  applicator  was  designed  and 
fabricated.  The  applicator  was  made  of  polystyrene  with  loading  ports  for  nine  seeds.  The 
central  portion  of  the  applicator  was  open  and  has  a  dimension  large  enough  for  tumor  to 
grow.  Equal  source  strength  was  assigned  to  all  nine  seeds  in  order  to  minimize  the 
possible  confusion  of  handling  variable  source  strengths.  The  seeding  configuration  was 
optimized  to  produce  an,  as  uniform  as  possible,  dose  distribution  to  the  central  portion  of 
the  applicator  and  to  be  usable  for  both  125I  and  103Pd  seeds.  Details  are  described  in 
Appendix  VI. 


Experiments  were  conducted  using  125I  seeds  with  an  initial  dose  rate  of  8  cGy/hr. 
Later  experiments  were  conducted  at  16  cGy/hr  using  1-125  alone,  103Pd  alone  and  a 
mixture  of  12T  and  103Pd  seeds.  Typical  results  are  shown  in  Figures  below.  Details  are 
described  in  a  manuscript  entitled  “Development  of  a  Rat  Solid  Tumor  Model  for 
Continuous  Low  Dose  Rate  Irradiation  Studies  using  125I  and  103Pd  Sources”  which  has 
been  prepared  for  submission  to  Int  J  of  Brachytherapy  Appendix  VI). 
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In  vivo  studies  using  an  in  vitro  assay 
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In  this  experiment,  BA1 112  tumor  cells  were  irradiated  in  vivo  to  graded  doses  from  2  to 
20  Gy.  Following  irradiation  at  low  dose  rate,  the  tumours  were  removed  and  their  colony 
formation  ability  was  measured  using  our  in  vitro  assay  techniques.  Figure  1 1  shows  a 
comparison  of  the  surviving  fraction  as  a  function  of  dose  calculated  for  the  BA1 112 
cells  to  that  measured  in  an  in  vivo/in  vitro  experiment  using  CLDRI 125I  irradiation.  The 
solid  line  represent  the  calculated  survival  curve  using  the  average  initial  dose  rates  used 
in  the  experiments.  The  open  circle  represent  calculation  using  the  actual  initial  dose  rate 
for  each  experiment.  The  parameters  of  a,  P,  repair  half-time,  and  tumor  doubling  time 
were  determined  directly  from  the  measurements  performed  on  the  BA1 1 12  cells  (details 
are  provided  separately). 
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Nude  Mice  Studies 


One  of  the  specific  aims  in  the  proposal  was  to  test  the  therapeutic  effectiveness 
of  Pd,  125I  and  a  Pd/I  mixture  in  human  prostate  carcinoma  xenografts  in  nude  mice, 
using  a  slow  growing  and  a  fast  growing  carcinoma.  This  was  to  be  followed  depending 
upon  the  successful  completion  of  the  BA1 12  in  rats  studies  and  to  be  initiated  near  the 
end  of  the  project.  Some  experimental  design  studies  were  performed  to  investigate  the 
feasibility  of  using  the  applicator.  Several  different  designs  were  considered  from  a 
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dosimetric  point  of  view.  Because  of  the  smaller  size  of  mice  and  the  limitations  in 
reducing  the  applicator  size  and  weight,  we  were  unable  to  reach  a  practical  solution. 
From  these  initial  studies,  we  have  concluded  that  xenografts  studies  on  nude  mice  would 
not  be  feasible  at  this  stage.  A  different  animal  model  and  brachytherapy  applicator 
would  be  needed  to  investigate  this  issue  further. 

KEY  RESEARCH  ACCOMPLISHMENTS 

•  A  theoretical  model  based  on  incomplete  repair  during  CLDRI  has  been 
developed  for  addressing  the  questions  raised  in  the  project.  The  BED  for 
implants  with  a  mixture  of  two  radionuclides  has  been  derived  as  an  analytical 
expression.  A  manuscript  describing  this  work  has  been  published  in  International 
Journal  of  Radiation  Oncology  (Appendix  I), 

•  Cell  survival  curves  for  both  125I  and  103Pd  were  measured  using  monolayers  of 
Chinese  hamster  cells  in  a  petri  dish  irradiated  at  low  dose  rates  using  12T  and 
103Pd  sources.  The  dose  sparing  effect  of  7  cGy/hr  relative  to  12  cGy/hr  can  be 
expressed  by  dose  modifying  factors  of  2  ±0.6  and  1.5  ±0.5  for  103Pd  and  125I, 
respectively.  The  RBEs  of  103Pd  relative  to  125I  were  1.2  ±  0.4  and  2.0  ±  0.5  for  7 
and  12  cGy/hr,  respectively.  In  our  system,  the  RBE  of  103Pd  at  19.7  cGy/hr 
relative  to  125I  at  7.72  cGy/hr  is  estimated  to  be  3±1.  A  manuscript  describing  this 
work  is  in  preparation. 

•  Cell  survival  curves  for  103Pd  were  measured  using  monolayers  of  B A1 1 1 2  cells 
in  a  petri  dish  irradiated  at  low  dose  rates  using  103Pd  sources.  An  orthovoltage  x- 
ray  machine  was  adapted  to  produce  nearly  monoenergetic  21  keV  photons, 
which  simulates  103Pd  photon  energies.  Using  this  x-ray  beam,  cell  survival 
curves  for  the  BA1 1 12  cells  in  a  petri  dish  irradiated  at  an  acute  dose  rate  were 
also  measured.  We  have  successfully  tested  our  theoretical  model  for  predicting 
the  CLDRI  survival  curves  in  vitro  for  103Pd  sources  from  the  acute  dose  rate 
exposure  data.  A  manuscript  entitled  “Dose  rate  dependence  of  the  relative 
biological  effectiveness  of  103Pd  for  Continuous  Low  Dose  Rate  Irradiation  of 
BA1 1 12  Rhabdomyosarcoma  Cells  in  vitro  relative  to  Acute  Exposures”  has  been 
prepared  for  submission  to  Int,  J.  Radiat,  Oncology.  (Appendix  II). 

•  Cell  survival  curves  for  125I  were  measured  using  of  BA1 1 12  cells  irradiated  in 
vivo  at  low  dose  rate  of  8  cGy/hr  using  the  afterloading  rat  applicator  with  125I 
sources.  An  orthovoltage  x-ray  machine  was  adapted  to  produce  nearly 
monoenergetic  28  keV  photons,  which  simulates  125I  photon  energies.  Using  this 
x-ray  beam,  cell  survival  curves  for  the  BA1 1 12  cells  in  a  petri  dish  irradiated  at 
an  acute  dose  rate  were  also  measured.  We  are  now  attempting  to  test  our 
theoretical  model  for  predicting  the  CLDRI  survival  curves  in  vivo  for  125I 
sources  from  the  acute  dose  rate  exposure  data. 

•  In  order  to  produce  a  consistent  dose  distribution  to  irradiate  the  tumors 
transplanted  to  different  animals  and  to  minimize  the  radiation  exposure  to 
personnel  handling  the  radioactive  seeds,  an  afterloading  seed  applicator  has  been 
designed.  Several  animals  have  been  treated  with  CLDRI  using  these  applicators 
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*  * 

containing  ~  I  seeds  at  9  cGy/hr.  Tumor  growth  was  significantly  slowed  by 
CLDRI,  from  a  tumor  doubling  time  of  2.7  days  in  controls  to  13  days  in  the 
treated  animals.  However,  no  tumor  cures  were  observed  at  8  cGy/hr.  Further 
experiments  were  conducted  at  16  cGy/hr  for  125I  alone,  103Pd  alone  and  a  50:50 
mixture  of  125I  and  l03Pd  seeds.  Tumor  cures  were  observed  at  the  higher  dose  rate 
in  all  cases.  A  manuscript  entitled  “Development  of  a  Rat  Solid  Tumor  Model  for 
Continuous  Low  Dose  Rate  Irradiation  Studies  using  125I  and  103Pd  Sources”  has 
been  prepared  for  submission  to  Int  J  of  Brachytherapy  (Appendix  VI). 


REPORTABLE  OUTCOMES 

A  manuscript  entitled  “Biologically  effective  dose  (BED)  for  interstitial  seed 
implants  containing  a  mixture  of  radio-nuclides  with  different  half  lives”  by  Zhe  Chen, 
Ph.D.  and  Ravinder  Nath,  Ph.D.  has  been  published  in  the  International  Journal  of 
Radiation  Oncology,  The  manuscript  is  attached  as  an  Appendix  I. 

A  manuscript  entitled  “Development  of  a  Rat  Solid  Tumor  Model  for  Continuous 
Low  Dose  Rate  Irradiation  Studies  using  125I  and  103Pd  Sources”  has  been  prepared  for 
submission  to  Int  J  of  Brachytherapy  (Appendix  VI). 

A  manuscript  entitled  “Dose  rate  dependence  of  the  relative  biological 
effectiveness  of  103Pd  for  Continuous  Low  Dose  Rate  Irradiation  of  BA1 1 12 
Rhabdomyosarcoma  Cells  in  vitro  relative  to  Acute  Exposures”  has  been  prepared  for 
submission  to  Int.  J.  Radiat.  Oncology.  (Appendix  II). 


CONCLUSIONS 

We  have  made  considerable  progress  towards  the  specific  aims  of  the  project. 
Theoretical  model  for  continuous  low  dose  rate  irradiation  using  a  mixture  of 
radionuclides  has  been  developed.  Experiments  have  been  performed  using  BA1112 
tumor  cells  and  Chinese  Hamster  cells  growing  in  vitro  and  BA1 1 12  cells  growing  in 
vivo  as  solid  tumors  in  WAG/rij  rats.  Radiobiology  parameters  for  these  cells  have  been 
determined  and  used  in  the  theoretical  radiobiology  model  to  improve  our  understanding 
of  the  experimental  observations.  We  have  designed  and  fabricated  applicators  for  in  vivo 
irradiations  as  well  as  developed  the  animal  care  procedures.  We  have  performed  in  vivo 
experiments  for  tumor  growth  studies  using  the  BA1 1 12  rat  model  with  125I ,  l03Pd  and  a 
50:50  mixture  of  the  two  seeds. 
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PHYSICS  CONTRIBUTION 


BIOLOGICALLY  EFFECTIVE  DOSE  (BED)  FOR  INTERSTITIAL  SEED 
IMPLANTS  CONTAINING  A  MIXTURE  OF  RADIONUCLIDES  WITH 

DIFFERENT  HALF-LIVES 


Zhe  Chen,  Ph.D.,  and  Ravinder  Nath,  Ph.D. 

Department  of  Therapeutic  Radiology,  Yale  University  School  of  Medicine,  New  Haven,  CT 

Purpose:  To  develop  a  tool  for  evaluating  interstitial  seed  implants  that  contain  a  mixture  of  radionuclides  with 
different  half-lives  and  to  demonstrate  its  utility  by  examining  the  clinical  implications  of  prescribing  to  an 
isodose  surface  for  such  an  implant. 

Methods  and  Materials:  A  linear-quadratic  model  for  continuous  low  dose  rate  irradiation  was  developed  for 
permanent  implants  containing  a  mixture  of  radionuclides.  Using  a  generalized  equation  for  the  biologically 
effective  dose  (BED),  the  effects  of  cell  proliferation  and  sublethal  damage  repair  were  examined  systematically 
for  implants  containing  a  mixture  of  radionuclides.  A  head-and-neck  permanent  seed  implant  that  contained  a 
mixture  of  125I  and  103Pd  seeds  was  used  to  demonstrate  the  utility  of  the  generalized  BED. 

Results:  An  equation  of  BED  for  implants  containing  a  mixture  of  radionuclides  with  different  half-lives  was 
obtained.  In  such  an  Implant,  the  effective  cell  kill  was  shown  to  depend  strongly  on  the  relative  dose 
contributions  from  each  radionuclide  type;  dose  delivered  by  radionuclides  with  shorter  half-life  always  resulted 
in  more  cell  kill  for  any  given  sublethal  damage  repair  and  cell  proliferation  rates.  Application  of  the  BED 
formula  to  an  implant  containing  a  mixture  of  12SI  and  i03Pd  seeds  demonstrates  that  the  conventional  dose 
prescription  to  an  isodose  surface  is  not  unique  for  such  an  implant.  When  the  prescription  dose  was  based  on 
existing  clinical  experience  of  using  125I  seeds  alone,  mixing  103Pd  seeds  with  125I  seeds  would  increase  the  cell 
kill.  On  the  other  hand,  if  the  prescription  dose  were  based  on  existing  clinical  experience  of  using  103Pd  seeds 
alone,  mixing  l25I  seeds  with  103Pd  seeds  in  the  same  implant  would  create  radiobiologically  “cold”  spots  (Le.,  an 
increase  in  cell  survival)  at  locations  where  a  major  portion  of  the  prescription  dose  is  contributed  by  the  125I 
seeds.  For  fast-growing  tumors,  these  “cold”  spots  can  become  significant. 

Conclusions:  Total  dose  alone  is  no  longer  sufficient  for  a  complete  characterization  of  a  permanent  seed  implant 
containing  a  mixture  of  radionuclides  with  different  half-lives  due  to  the  presence  of  cell  proliferation  and 
sublethal  damage  repair  in  the  protracted  dose  delivery.  BED  provides  a  tool  for  evaluating  the  radiobiologic 
effects  of  mixing  different  type  of  radionuclides  In  the  same  implant.  When  radionuclides  of  different  half-lives 
are  mixed  in  a  permanent  implant,  using  the  dose  prescription  established  from  existing  clinical  experience  of 
implants  with  the  longer  half-life  radionuclide  would  help  to  avoid  radiobiologic  “cold”  spots.  ©  2003  Elsevier 
Science  Inc. 

Interstitial  implant.  Biologically  effective  dose,  Iodine-125,  Palladium-103,  Brachytherapy,  Radiobiology. 
INTRODUCTION 

Permanent  implantation  of  encapsulated  radioactive  seeds 
in  tumors  has  been  used  widely  as  a  primary  or  adjuvant 
therapy  for  treating  prostate  and  head-and-neck  cancers 
(1-5).  At  present,  seeds  that  contain  the  radionuclide  of  l25I 
or  103Pd  are  routinely  used  in  permanent  interstitial  implants 
(6).  125I  seed  has  been  introduced  since  the  late  1960s  to 
replace  198 Au  seeds  due  to  its  long  half-life  [59.4  days  (7)], 
which  is  convenient  for  storage,  and  its  low  photon  energy, 
which  is  easy  for  radiation  protection.  It  has  remained  a 
popular  choice  for  permanent  interstitial  implant  (8).  103Pd 
seed,  which  has  an  average  photon  energy  close  to  l25I  and 


a  shorter  decay  half-life  [16.991  days  (7)],  was  introduced 
for  clinical  implant  about  a  decade  ago.  It  is  generally 
considered  that  103Pd  seeds,  due  to  their  short  half-lives,  are 
more  effective  for  fast  growing  tumors,  whereas  125I  seeds 
are  better  for  slow  growing  tumors  (9).  The  decision  regard¬ 
ing  which  radionuclide  to  use  for  a  given  implant  has  been 
influenced  largely  by  the  historical  development  of  radio¬ 
active  seeds  and,  recently,  by  the  radiobiologic  consider¬ 
ations  for  each  radionuclide  (9).  Nonetheless,  seeds  contain¬ 
ing  different  radionuclide  types  have  never  been  reported  in 
the  same  clinical  implant  For  tumors  that  may  contain  both 
fast  and  slow  growdng  cells,  it  would  seem  desirable  to  use 
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Fig.  1.  Accumulative  dose  delivered  as  a  function  of  implant  time  for  125I  and  103Pd  implants.  Note  that  the  vertical  axis 
was  plotted  as  the  ratio  of  delivered  dose  to  the  dose  to  full  decay. 


a  mixture  of  short  and  long  half-life  radionuclides  (e.g., 
103Pd  and  125I)  in  the  same  tumor  volume.  Mixing  seeds  of 
different  half-lives  in  the  same  implant,  however,  requires 
new  tools  to  evaluate  the  clinical  impact  of  the  resulting 
implant.  The  main  goal  of  this  study  is  to  present  a  tool  for 
evaluating  this  type  of  implant  and  to  demonstrate  its  utility 
by  examining  the  appropriateness  of  keeping  isodose  pre¬ 
scription  in  such  an  implant. 

The  need  for  a  tool,  other  than  the  delivered  dose,  stems 
from  the  well-known  observation  in  radiobiology  that  the 
cell  survival  in  a  permanent  interstitial  implant  depends  not 
only  on  the  total  dose  but  also  on  the  temporal  pattern  of 
dose  delivery  (9-15).  For  implants  containing  a  single  type 
of  radionuclide,  the  temporal  pattern  of  dose  delivery  is 
characterized  by  a  simple  exponential  reduction  of  dose  rate 
as  a  function  of  time.  Figure  1  illustrates  the  dose  delivered 
as  a  function  of  implant  time  for  a  125I  and  a  103Pd  implant. 
The  time  required  to  deliver,  for  example,  80%  of  dose  to 
full  decay  differs  by  a  factor  of  more  than  3  (about  140  days 
for  125I  implants  and  about  40  days  for  103Pd  implants).  This 
difference  in  the  rate  of  dose  delivery  (caused  by  the  dif¬ 
ferent  decay  half-lives)  can  result  in  very  different  clinical 
responses  for  a  given  total  dose  when  the  surviving  cells  in 
the  irradiated  volume  are  continuously  proliferating  and  the 
sublethally  damaged  cells  can  be  repaired  during  the  pro¬ 
tracted  dose  delivery  (9,  11,  12).  Indeed,  this  consideration 
has  led  the  prescription  dose  for  103Pd  seed  implants  in 
monotherapy  for  prostate  cancer  to  125  Gy,  considerably 
less  than  the  145  Gy  normally  prescribed  for  125I  implants. 
When  seeds  of  different  half-lives  are  mixed  in  the  same 


implant,  the  temporal  dose  delivery  pattern  to  a  point  within 
the  implant  is  now  defined  by  the  decay  characteristics  of  all 
radionuclide  types  that  contribute  doses  to  it.  The  resulting 
temporal  pattern  becomes  complex  (can  be  nonexponen¬ 
tial),  and  its  form  varies  throughout  the  implant  volume  (as 
a  function  of  seed  locations).  To  properly  evaluate  such  an 
implant,  a  tool  that  can  capture  the  interplay  between  the 
complex  temporal  dose  delivery  pattern  and  the  underlying 
tissue  radiobiology  is  needed. 

The  concept  of  biologically  effective  dose  (BED)  has 
been  used  as  a  surrogate  of  the  biologic  responses  to  ion¬ 
izing  radiation  (13-15).  For  a  course  of  fractionated  radio¬ 
therapy  or  brachytherapy,  the  surviving  fraction  of  the  irra¬ 
diated  cell  (S)  is  directly  linked  to  BED  as 

S  =  e"aBED  (1) 

where  a  characterizes  the  radiosensitivity  of  cells.  The 
linear-quadratic  (LQ)  cell  survival  model  has  been  used 
by  many  investigators  to  calculation  the  BED  for  acute, 
fractionated,  and  protracted  dose  delivery  (9-15).  For 
protracted  irradiation  at  a  constant  dose  rate,  Thames  has 
extended  the  LQ  model  for  acute  irradiation  to  take  into 
account  the  kinetics  of  sublethal  damage  repair  during  the 
dose  delivery  (13).  Later,  Dale  derived  an  analytic  for¬ 
mula  of  BED  for  continuous  irradiation  with  exponen¬ 
tially  decreasing  dose  rate  as  encountered  in  brachyther¬ 
apy  with  short-lived  radionuclides  (11,  12).  Dale’s 
formulation,  fundamentally  equivalent  to  the  incomplete 
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repair  model  (13),  takes  into  account  both  the  cell  pro¬ 
liferation  and  sublethal  damage  repair  during  the  dose 
delivery  characterized  by  a  simple  exponential  function. 
Ling  has  used  the  model  to  study  the  relative  radiobio¬ 
logic  effectiveness  (RBE)  of  permanent  implants  using 
198 Au,  125I,  or  ia?Pd  seeds  (9)  and  to  examine  the  biologic 
effects  of  dose  heterogeneity  inherent  in  interstitial  im¬ 
plants  (16).  In  this  article,  we  generalize  Dale’s  BED 
formula  for  implants  containing  a  mixture  of  radionu¬ 
clides  with  different  half-lives.  The  potential  of  using  the 
generalized  BED  formula  as  a  tool  for  evaluating  im¬ 
plants  containing  a  mixture  of  radionuclides  is  examined. 
Its  usage  is  demonstrated  by  examining  the  clinical  im¬ 
plications  of  isodose  prescription  for  a  head-and-neck 
implant  containing  a  mixture  of  125I  and  103Pd  seeds. 

METHODS  AND  MATERIALS 

Implants  containing  a  mixture  of  radionuclides 

For  simplicity,  let  us  consider  an  implant  containing  two 
types  of  radionuclides  (e.g.,  125I  and  I03Pd).  Assume  that 
there  are  Nt  seeds  of  radionuclide  type  1  and  N2  seeds  of 
radionuclide  type  2.  In  such  an  implant,  the  instantaneous 
dose  rate  to  a  point  r  at  time  t  after  implantation  is  given  by 

2 

D(r,t)  =  2  D0i(r)e~x “  (2a) 

r=l 

where  D0/(r)  denotes  the  total  initial  dose  rate  produced  by 
the  seeds  of  radionuclide  type  L  It  depends  on  the  spatial 
locations  of  all  type  i  seeds.  It  is  clear  from  Eq.  2a  that  the 
overall  temporal  pattern  in  such  an  implant  depends  on  the 
value  of  D0i(r)u  which  can  be  variable  throughout  the 
implant  volume.  For  point-like  sources,  D0/(r)  is  approxi¬ 
mately  equal  to  (17) 


£«Cr)  ~  8i(\r,i  ~  rlMJi  (2b) 

*=ilri7  n 

In  the  above  equations,  the  index  i  labels  the  radionuclide 
type  and  l  labels  the  individual  seeds.  Su  denotes  the  air- 
kerma  strength  for  the  seed  l  of  radionuclide  type  i.  Ah  Af, 
(<j>m h  and  8i(r)  denote  the  dose  rate  constant,  the  radioac¬ 
tive  decay  constant,  the  anisotropy  constant,  and  the  radial 
dose  function,  respectively,  for  the  seeds  of  radionuclide 
type  i .  r0  denotes  the  reference  distance  (usually  1  cm)  at 
which  the  dose  rate  constant  was  determined.  In  this  work, 
the  delivered  dose,  isodose  distribution,  and  dose-volume 
histograms  used  in  conventional  implant  dose  evaluation 
were  computed  by  using  Eq.  2  with  appropriate  parameters 
for  125I  and  103Pd  seeds  (17). 

Linear-quadratic  model  for  CLDRI 

Dale’s  work  on  BED  for  implants  containing  a  single 
type  of  radionuclide  (11,  12)  is  generalized  in  this  paper  for 
implants  containing  a  mixture  of  radionuclides  of  different 


half-lives.  To  model  the  kinetics  of  sublethal  damage  repair. 
Dale  has  invoked  the  assumption  that  radiation-induced  cell 
inactivation  is  caused  by  the  damage  of  two  critical  targets 
in  a  cell  to  model  the  kinetics  of  sublethal  damage  repair. 
Under  this  assumption,  when  a  radiation  event  damages 
only  one  critical  target,  the  cell  is  considered  sublethally 
damaged,  which  is  repairable.  Cell  kill  occurs  only  when  the 
other  critical  target  is  damaged  before  the  existing  damage 
is  fully  repaired.  By  assuming  sublethal  damage  repairs 
exponentially  with  time,  i.e.,  if  a  sublethal  damage  was 
inflicted  at  time  then  the  probability  for  it  persisting  to 
time  t  is  the  rate  of  sublethal  damage  repair  can 

be  characterized  by  a  single  parameter,  repair  half-time 
Ln  2 

T\n  ~  - *  Average  repair  half-times  reported  for 


mammalian  normal  and  tumor  tissues  vary  from  0.5  to 
3  h  with  normal  tissues  usually  possessing  longer  repair 
half-time  (18).  For  “generic”  tumors,  a  repair  half-time  of 
1.5  h  has  often  been  used  in  model  calculation  (19).  The 
rate  of  cell  proliferation  was  also  modeled  by  a  single 
parameter,  the  tumor  potential  doubling  time  Tp .  This 
model  assumes  that  the  cells  in  a  target  volume  all 
proliferate  at  the  same  rate  and  cell  loss  from  the  tumor 
volume  is  negligible.  The  tumor  potential  doubling  time 
is  tumor-tvpe  dependent  and  may  vary  from  patient  to 
patient  even  among  the  same  tumor  type.  For  squamous 
cell  head-and-neck  cancer,  a  “typical”  Tp  of  5  days  has 
been  quoted  in  literature  (19).  For  prostate  carcinoma,  Tp 
values  ranging  from  10  to  60  days  have  been  reported  in 
the  literature  (20).  Using  these  simple  models  for  suble¬ 
thal  damage  repair  and  cell  proliferation.  Dale  derived  an 
analytic  expression  of  BED  for  interstitial  implants  con¬ 
taining  a  single  type  of  radionuclide  as  follows. 


,  x  0.6937W 

BED  =  D(Teff)  RE - (3) 

where  RE  is  given  by 

\ 

s£=i+2(fe  <*•> 

with 

D\  f  1  1  1 

7  =  12A  (1  "  ~  AT H  (1  “ 

(4b) 

In  the  above  equations,  Z>0  denotes  the  initial  dose  rate, 
D0 

D{Teffi  =  —  (1  —  e  kTeff)  is  the  dose  delivered  up  to  Teffh 

and  a  and  jB  are  coefficients  of  the  LQ  model.  The  effect  of 
sublethal  damage  repair  is  captured  by  the  factor  RE,  and 
the  effect  of  cell  proliferation  is  characterized  by  the  second 
term  on  the  right-hand  side  of  Eq.  3.  This  analytic  form  of 
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BED  has  enabled  a  systematic  examination  of  the  interplay 
between  the  radioactive  half-life  and  the  underlying  suble- 
thal  damage  repair  and  cell  proliferation  rates.  Many  inves¬ 
tigators  have  discussed  the  use  of  this  BED  for  implants 
containing  a  single  type  of  radionuclide  (9,  11,  12,  16,  18, 
19, 21).  In  this  paper,  a  generalized  form  of  Eq.  4  is  derived 
for  implants  containing  a  mixture  of  any  number  of  radio¬ 
nuclides  with  different  half-lives. 


RESULTS 

BED  for  implants  containing  a  mixture  of  radionuclides 
with  different  half-lives 

Following  Dale’s  work,  a  generalized  BED  for  implants 
containing  a  mixture  of  radionuclides  with  different  half- 
lives  was  obtained  (see  Appendix  for  details  of  derivation). 
For  an  implant  mixed  with  two  different  half-lives,  the  y  in 
Eq.  4  is  replaced  by 

+  ^  {i (1  -  -  vb ( <>  - 

AojAb  f_J_ 

/*-  a21a,  +  a2u  ’ 

<s> 

where  D0l  and  D02  denote  the  initial  dose  rates,  at  the  point  of 
calculation,  produced  by  seeds  containing  radionuclide  type  1 
and  2,  respectively;  and 

D(Teff)  =  yji  (i—  C-A,r*)  +  ^  (1  - 


Equation  5  is  one  of  the  main  results  of  this  work.  Its 
general  property  and  application  to  a  clinical  implant  are 
discussed  in  the  remainder  of  this  section. 

BED  as  a  function  of  sublethal  damage  repair 
The  effects  of  sublethal  damage  repair  for  permanent 
implants  containing  seeds  of  different  half-lives  can  be 
examined  under  the  condition  of  Tp  =  <»,  i.e.,  no  cell 


proliferation.  In  this  case,  the  Teff  is  infinite  and  Eq.  5 
becomes 

j3 

RE  =  1  +  -  X 
a 

&02  A)lA)l(^l  +  A2  +  2fl) 

Ai(At  +  p)  A2(A2  +  /x)  (A1  +  A2)(A1  +  fji)(\2  +  p) 

A)2 

A2 

(6) 

for  an  implant  containing  two  different  half-lives.  BED  is 
simply  the  product  of  the  total  dose  and  RE.  Figure  2  plots 
the  RE  as  a  function  of  repair  half-time  for  an  implant 
containing  a  mixture  of  125I  and  103Pd  seeds  for  a  fixed  total 
dose.  Two  general  properties  are  conveyed  by  this  figure: 
First,  at  a  given  repair  half-time,  the  RE,  or  BED,  is  a 
function  of  relative  dose  contributions  from  the  two  types  of 
seeds.  By  varying  the  relative  dose  contribution  from  125I 
and  103Pd  seeds,  the  BED  for  the  mixed-seed  implant  can  be 
varied  continuously  from  that  of  using  125I  alone  to  that  of 
using  103Pd  alone.  Therefore,  the  BED  in  such  an  implant 
will  be  highly  dependent  on  the  locations  of  the  implanted 
seeds.  Second,  the  BED  is  larger  when  the  total  dose  is 
delivered  solely  by  the  radionuclide  with  shorter  half-life 
(103Pd  in  this  case).  This  point  is  further  illustrated  in  Fig.  3 
in  which  the  RE  associated  with  total  dose  delivered  by 
either  125I  or  103Pd  alone  is  plotted  from  7^  =  0  to  = 
00.  If  sublethal  damage  can  be  repaired  instantaneously  (i.e., 
7*1/2  “  0)  or  is  nonrepairable  (i.e.,  7^  =  00),  no  difference 
in  cell  kill  exists  between  different  radionuclides.  In  gen¬ 
eral,  implants  with  radionuclide  of  shorter  half-life  exhibit 
more  effective  cell  kill  for  the  same  delivered  dose.  This 
advantage  of  using  shorter  half-life  radionuclide  would  in¬ 
crease  initially  as  repair  half-time  increases  and  would 
diminish  eventually  when  sublethal  damage  becomes  irrep¬ 
arable. 

BED  as  a  function  of  cell  proliferation 
The  presence  of  cell  proliferation  would  always  reduce 
the  magnitude  of  BED,  because  part  of  the  delivered  dose 
has  to  counter  the  repopulated  cells.  The  amount  of  reduc¬ 
tion,  0.693t/(a7/?),  depends  on  both  the  proliferation  rate 
[0.693 f(ctTp)]  and  the  implant  duration  (/).  If  the  prolifera¬ 
tion  rate  remains  the  same  during  the  dose  delivery,  there 
would  exist  a  time  at  which  the  cell  inactivation  rate  in¬ 
duced  by  the  decaying  dose  rate  would  become  less  than  the 
proliferation  rate.  Beyond  this  time,  there  would  be  no  net 
cell  kill  from  the  remaining  dose  delivered.  The  time  at 
which  the  cell-kill  rate  equals  the  proliferation  rate  is  de¬ 
fined  as  the  effective  treatment  time  7^(12).  Figure  4  plots 
the  BED  calculated  at  Teff  as  a  function  of  the  potential 
tumor  doubling  time  and  the  relative  dose  contributions  in 
an  implant  with  mixed  125I  and  103Pd  seeds.  A  repair  half- 
life  of  1.5  h  and  an  a/p  of  10  were  used  in  the  calculation. 
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Fig,  2,  Relative  effectiveness,  RE,  calculated  as  a  function  of  repair  half-time  for  implants  using  a  mixture  of  125I  and 
103Pd  seeds.  The  BED  is  the  product  of  RE  and  the  total  delivered  dose. 


Total  delivered  dose  is  80  Gy.  When  the  relative  dose 
contribution  by  125I  and  103Pd  seeds  is  fixed,  the  BED  is 
shown  to  increase  with  Tp.  For  a  given  Tpt  BED  is  always 
larger  when  the  same  dose  is  delivered  by  the  radionuclide 
with  shorter  half-life  (103Pd  in  this  case).  The  difference  in 
BED  between  dose  delivered  by  103Pd  and  125I  alone  is  most 
significant  for  fast  growing  tumors  and  becomes  less  sig¬ 
nificant  for  slower  growing  tumors. 


Fig.  3.  Relative  effectiveness  (RE)  over  the  entire  range  of  possible 
repair  half-time  for  implants  using  125I  or  103Pd  alone. 


An  application  of  the  generalized  BED :  Inadequacy  of 
isodose  prescription  in  mixed  seed  implant 

The  basic  issue  regarding  the  isodose  prescription  for  an 
implant  containing  seeds  of  different  half-lives  is  as  follows. 
For  a  permanent  implant  using  only  a  single  type  of  radio¬ 
nuclide,  the  total  dose,  TD,  at  any  given  point  is  related  to 
the  initial  dose  rate,  D0, 

TD=  1.447W>0  (7) 

For  this  type  of  implant,  a  prescription  to  an  isodose  surface 
is  equivalent  to  a  prescription  to  an  isodose  rate  surface.  In 
other  word,  the  initial  dose  rate  is  fixed  once  a  prescription 
of  total  dose  is  established.  The  temporal  dose  delivery 
pattern  in  such  an  implant  follows  a  simple  exponential 
function  and  is  the  same  throughout  the  implant  volume. 
For  an  implant  containing  a  mixture  of  two  radionuclides, 
however,  Eq.  7  becomes 

TD  =  1 ,447f2D0(1)  +  1 .447 (8) 

The  total  dose  is  now  dependent  on  the  initial  dose  rates, 
Dq}  and  D(q\  produced  by  the  two  types  of  radionuclide, 
respectively.  According  to  Eq.  8,  a  prescription  to  an  isod¬ 
ose  surface  can  now  be  fulfilled  by  many  different  combi¬ 
nations  of  Dq}  and  IJ(q\  Therefore  the  prescription  to  an 
isodose  surface  is  not  unique  with  respect  to  the  initial  dose 
rates  produced  by  the  two  types  of  radionuclides.  Further- 
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Fig.  4.  BED  for  an  80  Gy  implant  using  a  mixture  of  125I  and  103Pd  seeds  plotted  as  a  function  of  tumor  potential 
doubling  time  Tp  and  the  fraction  of  total  dose  delivered  from  the  125I  seeds.  A  sublethal  damage  repair  half-life  of  1.5 
h  and  a  a/j8  of  10  for  tumor  were  used  in  the  calculation. 


more,  the  overall  temporal  pattern  of  dose  delivery  is  now 
variable  throughout  the  implant  volume,  determined  by  the 
sum  of  two  exponential  functions  with  variable  relative 
weightings. 

The  clinical  implications  of  prescribing  to  an  isodose 
surface  for  such  an  implant  can  be  examined  using  the 
generalized  BED  formula.  Let  us  take  a  clinical  implant 
performed  on  a  45-year-old  male  patient  who  presented 
with  cancer  of  the  larynx  as  an  example.  The  permanent 
seed  implant  to  the  postpharyngeal  wall  used  10  125I  seeds 
and  23  103Pd  seeds.  Initial  source  strength  was  0.64  U  for 
125I  seeds  and  0.88  U  for  103Pd  seeds.  Figure  5  shows  the 
isodose  lines  (to  full  decay),  in  the  x-y  plane  through  the 
center  of  the  implant,  contributed  by  the  125I  seeds  alone 
(top  panel),  the  103Pd  seeds  alone  (middle  panel),  and  by 
both  125I  and  103Pd  seeds  (bottom  panel).  The  80  Gy  isodose 
line  was  chosen  as  the  prescription  dose  for  the  implant. 
Figure  6  plots  the  BED  calculated  at  19  different  spatial 
locations  on  the  isodose  prescription  line  for  several  poten¬ 
tial  doubling  times.  The  BED,  dependent  on  the  relative 
dose  contribution  by  the  103Pd  and  125I  seeds,  is  different  at 
different  locations  on  the  same  isodose  line.  It  demonstrates 
that  the  isodose  prescription  is  inadequate  for  implants 
containing  seeds  of  different  half-lives.  Clinically,  if  the 
dose  prescription  were  established  from  existing  implant 
experience  that  uses  125I  seeds  only,  mixing  103Pd  seeds 
with  125I  seeds  would  increase  the  effective  cell  kill  for  the 
same  prescribed  dose.  On  the  other  hand,  if  the  dose  pre¬ 
scription  were  established  from  implant  experience  of  using 
103Pd  seeds  alone,  mixing  125I  seeds  with  103Pd  seeds  in  the 


same  implant  would  create  radiobiologically  “cold”  spots 
(an  increase  in  cell  survival  from  the  existing  clinical  ex¬ 
perience)  at  locations  where  a  major  portion  of  the  dose  is 
contributed  by  the  125I  seeds.  These  “cold”  spots  can  be¬ 
come  significant  for  fast  growing  tumors.  The  magnitude 
and  spatial  location  of  the  “cold”  spots  can  vary  throughout 
the  implant  volume. 

DISCUSSION 

As  shown  in  the  previous  sections,  new  tools  are  needed 
for  a  meaningful  implant  evaluation  when  radionuclides  of 
different  half-life  are  used  in  the  same  implant.  The  gener¬ 
alized  BED  formula  derived  from  Dale’s  work  (11,  12) 
provides  a  tool  to  examine  the  interplay  of  the  complex 
temporal  dose  delivery  pattern  in  such  an  implant  with  the 
underlying  radiobiologic  processes.  The  calculation  per¬ 
formed  in  this  work  suggests  that  the  BED  is  always  larger 
when  the  dose  is  delivered  by  103Pd  seeds  compared  with 
that  by  125I  seeds  for  any  given  sublethal  damage  repair  and 
cell  proliferation  rates  as  long  as  the  total  dose  is  the  same. 
Because  the  RBE  arising  from  different  linear  energy  trans¬ 
fer  (22),  secondary  electrons  generated  by  125I  and  103Pd 
photons,  is  not  dealt  with  in  this  paper,  the  above  conclusion 
is  due  purely  to  the  difference  in  decay  half-lives.  Because 
the  initial  dose  rate  is  always  greater  when  using  a  radio¬ 
nuclide  with  shorter  half-life,  the  conclusion  drawn  from  the 
BED  calculation  is  consistent  with  the  dose  rate  effect  (in 
the  absence  of  cell  cycle  effect)  observed  over  the  years  for 
protracted  dose  delivery  (10,  23),  albeit  the  dose  rate  is  now 
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Fig.  5.  Dose  distribution  on  the  x-y  plane  for  a  head-and-neck 
implant  using  l25I  and  103 Pd  seeds:  125I  seeds  alone  (top),  l03Pd 
seeds  alone  (middle),  and  from  the  actual  implant  with  mixed 
seeds  (bottom).  The  80  Gy  isodose  line  on  the  bottom  panel  was 
chosen  as  the  clinical  prescription  isodose. 


time-dependent  in  the  permanent  implants.  Parenthetically, 
in  vitro  measurements  have  shown  that  the  RBE  for  103Pd  is 
greater  than  that  for  ,25I  (22).  If  the  RBE  of  103Pd  and  125I 


were  taken  into  account  in  the  model  calculation,  the  ob¬ 
servations  made  in  this  article  would  have  been  further 
reinforced. 

Although  the  general  conclusions  drawn  from  the  model 
calculation  conform  to  common  intuition,  it  should  he  em¬ 
phasized  that  the  quantitative  numbers  of  BED  should  al¬ 
ways  be  viewed  in  the  context  of  the  limitations  inherent  in 
the  models.  Many  other  factors  that  affect  the  radiobiologic 
responses  of  human  tissues,  such  as  the  cell  cycle  effects 
(24),  radiation-induced  apoptosis  (25,  26),  and  tissue  archi¬ 
tecture,  were  not  considered  in  the  current  model.  Even 
then,  for  the  calculated  BED  to  be  clinically  meaningful, 
one  would  have  to  know  the  patient-specific  tumor  potential 
doubling  time,  sublethal  damage  repair  half-time,  and  other 
required  model  parameters.  Given  the  simplistic  nature  of 
the  biologic  model  and  the  lack  of  reliable  means  of  deter¬ 
mining  the  patient-specific  model  parameters  at  present,  it 
would  be  inappropriate  to  treat  the  calculated  BED  as  a 
quantitative  predictor  of  the  clinical  outcome  for  a  patient 
implant.  Taking  the  current  model  as  an  example,  the  rela¬ 
tive  uncertainty  in  the  calculated  BED  is  determined  by  the 
relative  uncertainty  of  each  model  parameter: 


ABED 

BED 


Aa  A(a/p)  ATt 

Air+B~^w+A 


-  +  C— — 


(9) 


where  the  coefficients  A,  /?,  and  C  are  given  by 


o.6937yrp 

D(Teff)RE-  0.693T^Tp’ 

(10a) 

D(Tejf){\  -  RE) 
D(TJRE  -  0.693T^Tp  ’ 

(10b) 

d  RE 

°<r->  a* 

C  D(Tejr)RE-0.693Tejj/T/ 

(10c) 

These  coefficients  determine  the  sensitivity  of  the  calculated 
BED  to  the  relative  uncertainties  that  may  incur  in  deter¬ 
mining  each  model  parameter.  Because  A,  By  and  C  are 
functions  of  the  model  parameters  {a,  a/p,  Tp ,  /*,},  the 
influence  of  the  relative  uncertainty  of  a  given  model  pa¬ 
rameter  on  the  calculated  BED,  in  fact,  depends  on  the 
characteristics  of  the  tumor  being  studied.  For  example,  the 
relative  uncertainty  associated  with  the  potential  tumor  dou¬ 
bling  time  has  a  negligible  effect  on  slowly  growing  tumors 
(i.e.,  A  — >  0,  for  Tp  — >  a>)  but  can  be  appreciable  on  fast 
growing  tumors.  Similarly,  relative  uncertainty  incurred  in 
determining  a/p  would  have  less  effect  on  tumors  with 
larger  a/p  than  on  tumors  with  small  a/p.  For  the  repair 
time  constant,  however,  BED  is  insensitive  to  Ajut/jui  when 
repair  is  extremely  slow  or  fast  and  is  most  sensitive  at 
intermediate  repair  kinetics  (see  coefficient  C  and  Fig.  3). 
Therefore,  the  numerical  values  of  the  calculated  BED  can 
not  and  should  not  be  taken  as  a  quantitative  indicator  of  the 
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Spatial  location  on  the  80  Gy  isodose  line 

Fig.  6.  BED  at  19  spatial  locations  along  the  80  Gy  prescription  isodose  line  for  the  head-and-neck  implant  with  l25I 
and  103Pd  seeds  for  various  tumor  potential  doubling  times.  A  sublethal  damage  repair  half-life  of  1.5  h  and  a  a/|3  of 
10  for  tumor  were  used  in  the  calculation. 


clinical  outcome  of  a  given  implant  at  present.  Nonetheless, 
the  BED  can  be  useful  in  comparing  the  relative  merits  of 
different  physical  implant  strategies  in  a  given  implant  by 
using  a  consistent  set  of  radiobiologic  model  parameters.  It 
is  on  this  premise  that  we  regard  the  BED  formula  as  a 
useful  tool  for  implant  evaluation. 

Many  interesting  issues  arising  from  mixing  seeds  of 
different  half-lives  in  the  same  implant  can  be  examined 
using  the  generalized  BED.  For  example,  as  mentioned  in 
the  “Introduction,”  it  is  interesting  to  examine  whether  a  net 
gain  in  local  control  can  be  achieved  by  implanting  both  125I 
and  103Pd  seeds  in  tumors  that  may  contain  both  fast  and 
slow  growing  cells.  If  the  answer  were  affirmative,  there 
may  exist  an  optimal  mixture  for  a  known  proportion  of 
slow  and  fast  growing  cells.  Also,  the  cells  in  the  surround¬ 
ing  normal  structures  are  known  to  exhibit  very  slow  pro¬ 
liferation  rates.  Because  radiobiologic  “cold”  spots  are  de¬ 
sirable  in  the  normal  structures,  it  would  be  interesting  to 
find  out  if  the  spatial  location  of  the  seeds  of  different 
half-lives  can  be  so  arranged  such  that  the  cell  survival  in 
the  nearby  normal  structures  is  maximized,  without  com¬ 
promising  the  effective  cell  kill  in  the  tumor  volume.  De¬ 
tailed  discussion  of  these  and  other  issues  is,  unfortunately, 
outside  the  scope  of  this  paper,  and  they  will  be  dealt  with 
elsewhere. 

The  clinical  implant  used  for  examining  the  appropriate¬ 
ness  of  isodose  prescription  in  implants  containing  seeds  of 
different  half-lives  is  one  of  two  such  implants  ever  per¬ 


formed.  Because  the  dose  prescription  for  the  implant  was 
based  on  the  experience  of  using  125I  seeds  alone,  the 
effective  cell  kill  in  the  clinical  implant  is  expected  to  be 
greater  than  those  implants  using  125I  seeds  alone.  These 
two  implants,  however,  cannot  be  used  to  answer  a  more 
general  question,  i.e.,  whether  the  expected  increase  in  cell 
kill  will  translate  to  increased  local  control  when  compared 
with  existing  implant  experience  using  125I  alone.  Such  a 
question  can  only  be  answered  through  systematic  clinical 
trials  on  a  large  number  of  patients  with  known  clinical 
presentations. 

For  a  complete  evaluation  of  an  implant  containing  seeds 
of  different  half-lives,  BED  over  the  entire  volume  of  in¬ 
terest  should  be  examined  because  the  total  dose  and  initial 
dose  rate  are  inherently  inhomogeneous.  With  the  advance 
of  medical  imaging  modalities  and  image-based  treatment 
planning  (27),  we  hope  that  three-dimensional  tissue  struc¬ 
tures  can  be  identified  more  accurately  and  that  a  3D  dis¬ 
tribution  of  BED  can  be  calculated  with  confidence  (28). 
With  the  BED  distribution,  radiobiologically  significant 
“cold”  spots,  due  either  to  “cold”  spots  in  total  dose  or  to  the 
“cold”  spots  in  dose  rate,  can  be  identified  during  the 
planning  of  an  interstitial  seed  implant. 

CONCLUSIONS 

A  generalized  Dale  equation  for  BED  is  presented  for 
interstitial  implants  containing  a  mixture  of  radionuclides 
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with  different  half-lives.  The  generalized  BED  provides  a 
tool  for  evaluating  the  radiobiologic  effects  of  mixing  dif¬ 
ferent  types  of  radionuclides  in  the  same  implant.  Model 
calculation  performed  in  this  paper  suggests  that  adding  the 
103Pd  to  the  125I  implant  would  increase  the  effectiveness  of 


cell  kill,  whereas  the  opposite  is  not  true  if  the  dose  pre¬ 
scription  was  based  on  the  clinical  experience  established 
with  the  125I  implants.  It  is  hoped  that  this  work  will 
stimulate  further  research  interest  in  improving  the  radiobi¬ 
ologic  modeling. 
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APPENDIX 


Equation  5  can  easily  be  derived  following  the  work  of 
Dale  (11,  12)  for  implants  containing  a  single  type  of 
radionuclide.  The  reader  is  referred  to  Dale’s  original  arti¬ 
cles  (11,  12)  for  the  detailed  steps  and  the  rationale  under¬ 
lying  the  derivation.  The  main  derivation  steps  are  outlined 
below. 

Dale  has  invoked  the  assumption  that  a  lethal  radiation 
damage  is  caused  by  the  damage  of  two  critical  targets  in  a 
cell.  When  the  two  critical  targets  are  damaged  simulta¬ 


neously  by  a  radiation  event,  the  resulting  lethal  damage  is 
termed  type  A  damage.  When  a  radiation  event  damages 
only  one  critical  target,  the  cell  is  considered  subiethally 
damaged,  which  is  repairable.  In  the  latter  case,  a  lethal 
damage  results  when  the  second  critical  target  is  damaged 
before  the  existing  sublethal  damage  is  fully  repaired.  This 
type  of  lethal  damage  is  termed  type  B  damage.  Type  A 
damage  is  always  proportional  to  the  total  delivered  dose 
irrespective  of  dose  rate,  whereas  the  type  B  damage  is  dose 
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rate  dependent  because  the  independently  damaged  targets 
may  be  repaired  over  time.  The  overall  radiation-induced 
cell  inactivation  is  therefore  proportional  to  the  sum  of  type 
A  and  type  B  damages.  The  main  goal  in  deriving  the  BED 
is  to  determine  the  total  probability  of  type  A  and  type  B 
damages  over  the  time  period  of  a  given  dose  delivery. 

For  an  implant  with  a  mixture  of  radionuclides,  the  in¬ 
stantaneous  dose  rate  to  a  point  of  interest  at  time  t  after  the 
implant  is  given  by  Eq.  2  in  the  main  text  as 

./>(*)  =  5  (Al) 

1=1 

where  the  summation  is  over  different  types  of  radionu¬ 
clides  used  in  the  implant.  Using  this  expression  for  instan¬ 
taneous  dose  rate,  the  derivation  for  BED  follows  almost 
exactly  the  same  steps  outlined  in  the  appendix  of  Dale’s 
article  (11,  12). 

The  type  A  damage  is  proportional  to  the  dose  delivered. 
Therefore,  the  total  type  A  damage  accumulated  over  a 
period  of  time  T  after  seed  implantation  is  given  by 


In  the  above  expressions,  a  and  J3  are  the  linear  and  qua¬ 
dratic  coefficients,  respectively,  of  the  linear-quadratic  cell 
inactivation  model,  p  is  the  repair  time  constant  that  char¬ 
acterizes  the  rate  of  sublethal  damage  repair. 

With  Eqs.  A2  and  A3,  the  ratio  of  total  lethal  damage  to 
the  type  A  lethal  damage  alone,  termed  relative  effective¬ 
ness  by  Dale,  can  be  determined  as 


RE  =  1  + 


Total  type  -  B  damage 
Total  type  -  A  damage 


(A4) 

Equation  A4  is  the  main  result  of  this  derivation.  The  BED 
for  an  implant  with  a  mixture  of  radionuclides  is  therefore 
given  by 


Total  type  -  A  damage 


e~KiT) 


(A2) 


The  type  B  damage  is  determined  by  the  joint  probability  of 
two  targets  being  damaged  at  different  times  by  two  differ¬ 
ent  radiation  events.  By  assuming  the  sublethal  damage 
repairs  exponentially  with  time,  the  total  type  B  damage 
accumulated  over  the  same  period  of  time  T  after  seed 
implantation  can  be  determined. 


.  •  _  0.693 

BED  =  (1  -  e-^RE  -  Teff  (A5) 

where  Tp  is  the  tumor  potential  doubling  time.  Note  that  Tejf 
is  the  effective  treatment  time  at  which  the  tumor  cell- 
proliferation  rate  exceeds  the  tumor  cell-inactivation  rate 
caused  by  the  instantaneous  dose  rate  at  that  time. 

As  expected,  for  implants  with  a  single  type  of  radionu¬ 
clide,  i  =  j  =  1,  and  Eqs.  A4  and  A5  reduce  to  Dale’s 
original  equation 
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Abstract 


Relative  biological  effectiveness  (RBE)  of  continuous  low  dose  rate  irradiation  using 
Pd  sources  was  measured  relative  to  acute  exposures  using  an  x-ray  beam  which  was 
simulated  to  generate  a  nearly  monoenergetic  photon  energy  similar  to  that  of  103Pd.  The 
survival  curves  were  fitted  to  the  linear-quadratic  model  for  continuous  low  dose  rate 
irradiation.  A  profound  dose  rate  effect  is  observed  at  the  low  dose  rates  in  the  range  of 
6.8  to  14.4  cGy/hr  that  are  typical  of  permanent  interstitial  brachytherapy.  At  a  survival 
level  of  0.01,  the  RBE  of  CLDRI  at  6.8  and  14.4  cGy/hr  using  10  Pd  sources  relative  acute 
exposures  at  a  dose  rate  of  333  cGy/hr  from  a  simulated  x-ray  beam,  is  reduced  by  a 
factor  of  3  and  2,  respectively. 


Methods  and  Materials 

BA1112  cells 


Animals  were  implanted  with  transplanted  BA1 1 12  rhabdomyosarcomas  by  the 
inoculation  of  7500  tumor  cells,  suspended  in  0.05  ml  of  sterile  cell  culture  medium,  into 
a  subcutaneous  site  on  the  head.  The  tumors  were  allowed  to  grow  for  3  weeks,  to  an 
experimental  volume  of  approximately  100-200  mm3.  The  animals  chosen  for  the  in  vitro 
irradiations  of  the  BA1 112  sarcoma  cells  were  euthanized  by  anesthetic  overdose  and  the 
tumor  was  removed  using  aseptic  techniques.  The  tumor  was  chopped  into  a  fine  mash 
with  a  sterile  blade,  and  the  tumor  pieces  added  to  a  trypsin  solution  and  stirred  at  37°C 
for  15  minutes.  The  suspension  is  now  filtered  to  remove  pieces  of  contaminants  and 
centrifuged  at  400G  for  10  minutes,  then  re-suspended  in  10ml  of  DMEM  culture  media. 
This  single-cell  suspension  of  tumor  cells  will  be  counted,  and  assayed  for  viability  using 
the  same  colony  formation  assay  used  for  cells  in  cultures  and  seeded  into  experimental 
dishes  at  cell  concentrations  of  5xl03  to  1.5xl06.  The  dishes  are  put  into  a  37°C 
incubator  with  95%  air/5%  C02.  These  cells  are  allowed  to  attach  to  the  petri  dishes 
(60mm),  and  begin  to  divide,  at  least  48  hrs,  before  dishes  are  used  for  inv  vitro 
irradiations. 


In  vitro  irradiation  technique  for  acute  irradiation 

The  radiobiological  response  of  the  cells  under  acute  exposure  condition  was 
studied  using  a  simulated  x-ray  beams  with  average  energies  equivalent  to  that  emitted 
Pd  (20  -  22.7  keV  with  an  average  of  20.5  keV)  generated  by  an  orthovoltage  unit 
(Pantek,  CT).  The  simulated  beams  not  only  have  the  average  energies  similar  to  that 
given  by  the  radioactive  isotopes  but  also  have  a  narrow  photon  energy  spectrum.  The 
narrow  photon  energy  spectrum  was  achieved  by  optimizing  the  tube  voltage  (which 
determines  the  upper  limit  of  the  photon  energy)  and  by  added  filtration  (which  filters  out 
the  low-energy  Bremsstrahlung  photons),  following  the  work  of  Muench  et  al5. 

The  beam  quality  of  the  simulated  beam  was  monitored  while  the  tube  voltage 
and  the  added  filtration  were  optimized.  The  half-value-layer  of  the  simulated  beam  was 
measured  using  aluminum  filters  (with  thickness  of  0.1  to  1.0  mm)  from  Nuclear 
Associates  under  narrow  beam  condition.  A  narrow  circular  beam,  with  a  diameter  of  6 
cm  at  source  to  chamber  distance  (SCD)  of  50  cm,  was  produced  a  homemade  lead 
collimator  mounted  to  DXT  300’s  accessory  mount.  The  aluminum  sheets  added  to  the 
beam  for  the  HVL  measurement  were  taped  to  the  bottom  of  the  lead  collimator.  An  air- 
equivalent  Spokas  chamber  (Exradin,  Model  No:  Al  (0.5  ml,  AE  plastic))  was  used  to 
measure  the  ionization  at  the  SCD  of  50  cm  in  air.  Due  the  energy  dependence  of  the 
chamber  at  the  low  energies,  the  measured  ionization  were  converted  to  corresponding 
exposure  and  the  HVL  is  then  determined  from  the  relative  exposure  as  a  function  of 
aluminum  filter  thickness.  Figure  1  plots  the  relative  exposure  as  a  function  of  filter 
thickness  obtained  for  the  HVL  measurement.  The  measured  HVL  was  then  used  to 
determine  the  energy  of  an  equivalent  mono-energetic  photon  beam.  The  HVL  as  a 


function  of  mono-energetic  photon  beam  energy  for  aluminum  is  taken  from  Johns  and 
Cunningham6.  The  optimized  103Pd  x-ray  beam  has  a  HVL  of  0.84  mm  Aluminum  and 
equivalent  mono-energetic  photon  energy  of  20.7  keV  (average  energy  from  103Pd 
sources  is  20.9  keV).  Table  I  summarizes  the  setting  for  the  simulated  103Pd  beam  used 
for  acute  in  vitro  cell  survival  studies. 


Radiation  output  for  the  simulated  beam  was  determined  by  the  measured  average 
dose  to  tissue  culture  medium  in  the  petri  dish  using  Fricke  dosimetry.  It  was 
crosschecked  using  the  air-equivalent  Spokas  chamber  and  a  calibrated  parallel-plate  soft 
energy  chamber.  Since  the  Fricke  dosimetry  was  carried  out  using  the  same  petri  dish 
with  similar  amount  of  solution  as  used  for  in  vitro  irradiation,  its  measured  output  was 
used  in  computing  the  dose  for  this  experiment.  The  dose  rate  to  the  petri  dish  was 
determined  to  be  33.3  Gy/hr. 

In  vitro  irradiation  technique  for  CLDRI 

A  specially  designed  103Pd  source  irradiators  made  of  polystyrene  with  a  central 
hole  for  placing  the  petri  dishes  on  top  of  the  sources,  with  or  without  some  polystyrene 
spacers  between  the  sources  and  the  tissue  culture  dishes  was  built  in  our  lab  (Fig.  1).  The 
irradiator  consists  of  a  20.3  x  20.3  x  10.2  cm  polystyrene  phantom  with  a  10.2  cm 
diameter  centered  hole  containing  a  polystyrene  disk  loaded  with  103Pd  sources.  The 
irradiator  was  loaded  with  26  103Pd  sources  (model  200;  Theragenics  Corp.,  Norcross, 
Atlanta),  each  with  an  initial  activity  of  xx  mCi.  The  sources  are  arranged  in  a  circle  with 
a  diameter  of  6  cm,  and  two  sources  are  placed  next  to  each  other  at  the  center  of  the 
circle.  These  sources  were  contained  in  a  1.25  cm  thick,  10  cm  diameter  polystyrene  disc. 
Different  dose  rates  spanning  the  range  from  6  to  13  cGy/hr  were  obtained  by  varying  the 
thickness  of  the  polystyrene  spacers.  The  103Pd  irradiator  was  placed  in  a  37°C  water- 
jacketed  incubator  and  surrounded  in  lead  foil  of  1  mm  thickness  to  shield  the  photons. 

The  dose  to  cells  in  the  monolayers  was  determined  from  the  measured  average 
dose  to  tissue  culture  medium  in  the  petri  dish  using  Fricke  dosimetry,  with  a  calculated 
correction  for  interface  effects  due  to  photoelectric  effect  in  the  tissue  culture  dishes  []. 
The  uniformity  of  dose  rate  across  the  petri  dish  was  verified  by  three  independent 
means:  LiF  thermoluminescent  dosimetry,  film  dosimetry  and  dose  calculations  by  a 
computerized  treatment  planning  package  (Theraplan).  In  all  cases,  the  dose  uniformity 
across  the  petri  dish  was  better  than  ±5%. 

Cytotoxicity  studies 

All  experiments  were  performed  using  cell  monolayers  that  were  prepared  by 
plating  cells,  suspended  from  exponentially  growing  stock  cultures,  into  60-mm  diameter 
Falcon  tissue  culture  dishes.  The  cells  were  incubated  for  18  hr  prior  to  irradiation,  to 
allow  them  to  attach  and  to  progress  into  exponential  growth.  Once  the  cells  were  in 
exponential  growth,  the  growth  medium  was  removed  and  replaced  by  fresh  growth 
medium  just  before  the  beginning  of  the  irradiations,  which  lasted  1  to  60  hr  in  CLDRI. 


Cells  were  maintained  in  a  humidified  95%  air-5  %  C02  environment  at  37°C 
during  irradiations,  as  well  as  during  control  conditions.  After  irradiation,  the  cells  were 
washed  with  Hanks’  balanced  salt  solution,  trypsinized,  and  counted  using  a  Coulter 
counter,  equipped  with  a  Channelizer  to  allow  assessment  of  and  correction  for  any 
changes  in  cell  size.  Cells  were  plated  for  colony  formation  in  at  least  four  dishes  per 
data  point  and  were  allowed  to  grow  in  a  humidified  95%  air-5%  C02  environment  at 
37  C  for  10  days.  The  colonies  were  then  fixed,  stained,  and  counted.  Experiments  were 
planned  to  obtain  approximately  120  colonies  in  each  dish,  by  adjusting  the  number  of 
cells  plated  per  dish  appropriately.  Cytotoxicity  was  characterized  using  cell  survival 
fraction,  which  is  a  ratio  of  plating  efficiencies  of  the  treated  cells  relative  to  those  of 
unirradiated  cells  under  control  conditions.  In  these  experiments,  however,  the  numbers 
of  cells  in  cultures  treated  with  CLDRI  were  significantly  lower  than  the  numbers  of  cells 
in  control  cultures,  even  though  the  same  numbers  of  cells  were  used  to  initiate  the 
cultures,  because  of  reduced  rates  of  cell  division  in  the  irradiated  cultures  and  because 
some  cells  died  during  the  prolonged  irradiations.  The  surviving  fractions  were  therefore 
corrected  to  account  for  the  deficit  in  cell  number  in  the  irradiated  cultures,  using  the 
formula: 


S  = 


exptl 


control 


control 


where 


P expti  =  plating  efficiency  of  the  experimental  sample; 

Pcontroi  =  plating  efficiency  of  the  control  sample; 

Cexpti  =  number  of  cells  harvested  from  the  experimental  sample  at  the  end 
of  irradiation  time; 

Ccontroi  =  number  of  cells  harvested  from  the  control  sample  assayed  at  the 
end  of  irradiation  time  (without  irradiation). 


Plating  efficiencies  (colonies  formed/cells  plated)  for  untreated  cells  were 
approximately  90%  in  the  experiments  reported  here. 


RESULTS  AND  DISCUSSIONS 


The  cell  loss  in  irradiated  cultures,  described  in  the  materials  and  methods  section, 
was  calculated  as  (1  -  Cexpti/CControi)  for  all  dose  rates  studied  here  (Figure  3).  These  data 
were  fitted  with  straight  lines  with  no  intercepts.  The  slopes  of  these  lines,  i.e.  the  rates 
of  cell  loss  during  irradiation,  increased  monotonically  with  increasing  dose  rate  for  both 
cell  lines  and  both  radioisotopes.  The  slopes  were  similar  for  the  two  cell  lines.  This 
confirms  that  the  cell  loss  is  basically  a  dose  rate  effect  primarily  related  to  the  length  of 
irradiation  period. 


Figure  4  plots  the  cell  survival  curves  for  the  BA1 1 12  cells  irradiated  by  103Pd  at 
different  dose  rates.  The  lines  are  fit  to  the  measured  data  using  the  linear-quadratic  (LQ) 
survival  model.  At  low  dose  rate,  the  survival  data  were  fitted  only  to  the  linear  term  of 
the  LQ  model  because  the  quadratic  term  did  not  significantly  improve  the  goodness  of 
the  fit.  If  the  fit  to  the  low  dose  rate  survival  data  can  be  extrapolated  to  higher  dose  (not 
guaranteed),  the  relative  biological  effectiveness  (RBE)  of  the  CLDRI  can  be  compared 
to  the  acute  irradiation  using  the  simulated  103Pd  beam.  At  the  survival  level  of  0.01,  the 
RBE  of  CLDRI  at  dose  rates  of  14.4  cGy/hr  and  6.8  cGy/hr  was  found  to  be  0,49  and 
0.31,  respectively.  The  relative  effectiveness  as  compared  to  acute  irradiation  was 
reduced  by  a  factor  of  2  and  3  for  the  CLDRI  dose  rates  of  14.4  cGy/hr  and  6.8  cGy/hr, 
respectively. 


(Model  calculation  will  be  added  later) 


Conclusion 

Relative  biological  effectiveness  of  continuous  low  dose  rate  irradiation  using 
Pd  sources  has  been  measured  relative  to  acute  exposures  using  an  x-ray  beam  which 
was  simulated  to  generate  a  nearly  monoenergetic  photon  energy  similar  to  that  of  103Pd. 
A  profound  dose  rate  effect  is  observed  at  the  low  dose  rates  in  the  range  of  6.8  to  14.4 
cGy/hr  that  are  typical  of  permanent  interstitial  brachytherapy.  At  a  survival  level  of  0.01, 
the  RBE  of  CLDRI  at  6.8  and  14.4  cGy/hr  using  103Pd  sources  is  reduced  by  a  factor  of  3 
and  2,  respectively,  relative  acute  exposures  at  a  dose  rate  of  333  cGy/hr  from  a 
simulated  x-ray  beam. 


Table  I.  Beam  characteristics  of  simulated  103Pd  x-ray  beams 


Beam 

Pd-103 

Equivalent 


kV 

mA 

Added  Filter 
(mm  AL) 

Beam  HVL 
{mm  AL) 

Energy 

Homogeneity 

<%) 

Equivalent 

Energy 

(keV) 

<E>  from 
isotope 
*  (keV) 

29 

25 

1.826 

0,84 

88.6 

20.7 

20.7 

Table  II RBE  of  CLDRI  using 

103Pd  sources 

Radiation  Source 

Dose  Rate 
(Gy/hr) 

RBE  relative  to 
simulated  x-rays 

Pd-103 

0.068 

0.49 

Pd-103 

0.144 

0.31 

Simulated  x-rays 


33.3 


1 


laitive  Exposure 


Aluminum  (mm) 


Figure  1 


Figure  2. 103Pd  irradiator 


Figure  3.  Cell  loss  curves 
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NIST2000 
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APPENDIX  III 

Radiobiological  parameters  for  the  CCL-16  cell  line 


1.  Tumor  growth 

a.  In  vitro  tumor  cell  growth 

Method: 

Results:  Figure  1  plots  the  in  vitro  tumor  cell  growth  as  a  function  of  time.  The  cell 
doubling  time  for  tumor  proliferation  is  estimated  by  fitting  the  cell  growth  to  an 
exponential  function.  The  was  found  to  be  13  hours  (0.55  days)  (file: 
CCL_Growth.xls) 

2.  Sublethal  damage  repair  time  constant 
a.  Split  dose  experiments 

Method:  Split  dose  experiment  was  performed  for  CCL-16  cell  line  using  the  250  kV  x- 
rays.  A  total  dose  of  8  Gy  was  given.  The  survival  fraction  was  measured  as  a  function  of 
time  interval  between  two  4  Gy  irradiation.  The  cells  were  kept  at  22°C  during  the 
waiting  interval  to  prevent  proliferation.  According  to  incomplete  repair  model  we  have, 

ln[S(Q/S(0])  =l 
ln[S(°o)/S(0)]  ^ 

where  S(t)  is  the  surviving  fraction  with  split  time  interval  t  and  (I  is  the  sub-lethal 
damage  repair  time  constant. 

Results:  Figure  2  plots  the  measured  survival  and  the  curved  fitted  to  the  above  equation. 
The  repair  half-time  was  found  to  be  13  minutes  (file:  ccl_splitdose.xls). 

3.  a  and  (3  for  the  linear-quadratic  cell  survival  model 

The  parameters  of  a  and  (3  for  the  CCL-16  cell  line  for  1251,  103Pd,  and  250  kV  photons 
were  determined  from  the  survival  curves  obtained  under  acute  exposure  condition  using 
simulated  125I  and  103Pd  x-ray  beam  on  an  orthovoltage  unit.  To  obtain  the  a  and  (3,  the 
acute  survival  curves  for  the  simulated  1251,  103Pd,  and  the  250  kV  x-rays  were  fitted  to 
the  LQ  cell  survival  equation.  Two  fitting  procedures  were  used.  In  the  first  procedure, 
the  cx/p  ratio  were  kept  the  same  for  all  three  radiation  beam  qualities.  The  a  for  each 
beam  quality  and  0(7 (3  were  treated  as  fitting  parameters  in  the  global  fit  (Fig.3A).  The 
second  approach  is  similar  to  the  first  procedure  except  the  (3  was  regarded  the  same  for 
all  three  beam  qualities  (Fig.3B).  The  fitted  a  and  (3  is  given  in  the  attached  table.  It 
should  be  noted  that  the  relationship  between  the  cell  survival  curve  of  250  kV  x-rays  and 
those  from  125I  and  103Pd  does  not  appear  consistent  with  the  general  trend  observed  on  the 
cell  kill  as  a  function  of  LET.  Investigation  on  this  observation  is  being  conducted.  Due 
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to  this  observed  inconsistency,  a  global  fit  using  only  the  125I  and  103Pd  beams  were  also 
performed  (Fig.  4A  with  fixed  a/ (3  and  Fig.4B  with  fixed  (3) .  The  fitted  parameters  are 
listed  in  a  separate  table. 


Radiobiological  Parameters  Determined  from  Acute  Exposure  Experiments 


Beam  quality 

Parameters  for  CCL-16 

a 

(Gy’1) 

P 

(Gy'2) 

a/p 

(Gy) 

Tumor 

Doubling  Time 
(day) 

Sub-lethal 
Damage  Repair 
Half-time  (hr) 

250  kV 

0.17 

0.037 

4.7 

0.55 

0.21 

I  equivalent 

0.12 

0.025 

4.7 

0.55 

0.21 

103Pd  equivalent 

0.16 

0.033 

4.7 

0.55 

0.21 

250  kV 

0.23 

0.031 

7.4 

0.55 

0.21 

125I  equivalent 

0.07 

0.031 

2.1 

0.55 

0.21 

103Pd  equivalent 

0.18 

0.031 

5.7 

0.55 

0.21 

Radiobiological  Parameters  Determined  from  Acute  Exposure  Experiments 

(using  only  12SI  and  103Pd  beams) 


Beam  quality 

Parameters  for  CCL-16 

a 

(Gy1) 

P 

(Gy'2) 

a/p  Tumor 

.  Doubling  Time 

(Gy)  «w> 

Sub-lethal 
Damage  Repair 
Half-time  (hr) 

125I  equivalent 

0.17 

0.021 

8.3 

0.55 

0.21 

103Pd  equivalent 

0.22 

0.027 

8.3 

0.55 

0.21 

125I  equivalent 

0.14 

0.024 

5.9 

0.55 

0.21 

103Pd  equivalent 

0.25 

0.024 

10.6 

0.55 

0.21 

Ln(S(t)/S(0))/Ln(S(inf)/S(0)) 
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AAPENDIX IV 

Radiobiological  parameters  of  the  BA1112  rhabdomyosarcomas 


1.  Growth  characteristics 

a.  In  vitro  cell  growth 

Methods:  In  vitro  tumor  cell  growth  rate  was  determined  by  examine  the  cell  growth, 
Ln{N(t)/N(0)},  as  a  function  of  time. 

Results:  Figure  1  plots  the  measured  Ln{N(t)/N(0)}  as  a  function  of  time  from  four 
experiments.  (Explanation  of  the  cell  number  reduction  at  large  times).  The  tumor  cell 
doubling  time  in  the  exponential  growth  phase  (between  the  elapsed  time  of  0  and  96 
hours)  was  found  to  be  20  hours  from  least  square  fit.  (file:  growth  curve  invitro  xc.xls) 

b.  In  vivo  tumor  growth 

Methods:  Animals  were  implanted  with  transplanted  B A1 1 12  rhabdomyosarcomas  by 
the  inoculation  of  7500  tumor  cells,  suspended  in  0.05  ml  of  sterile  cell  culture  medium, 
into  a  subcutaneous  site  on  the  head.  An  aseptic  technique  was  used  during  the 
inoculation.  The  rats  were  anesthetized  with  xylaxine/ketamine  (0.005  mg/g  and  0.075 
mg/g,  respectively)  given  IP,  and  a  small  area  on  the  rats’  head  was  shaved  for  the 
subcutaneous  injection  of  tumor  cells.  The  actual  inoculation  of  tumor  cells  take 
approximately  3  minutes  followed  by  the  administering  of  Yohimbine  (0.002  mg/g)  IP  to 
counteract  the  effects  of  the  xylazine.  The  animal  is  awake  and  alert  in  approximately 
20-30  minutes.  The  tumors  were  allowed  to  grow  tumor  volume  were  measured  as  a 
function  of  elapsed  time.  The  tumor  volume  is  calculated  as  an  ellipsoid  with  a,  b,  c 
measured  by  calipers. 

Results:  Figure  2A  plots  the  measured  tumor  volume  as  a  function  of  time.  The  tumor 
volume  growth  exhibited  two  growth  phases:  an  initial  phase  with  an  apparent  tumor 
volume  doubling  time  of  about  1.1  days  (or  26  hours  slightly  longer  than  the  in  vitro  cell 
doubling  time  of  20  hours)  and  a  late  “slow”  phase  with  an  apparent  volume  doubling 
time  of  2.7  -  2.8  days.  The  late  volume  growth  is  significantly  slower  than  the  tumor  cell 
doubling  time,  reflecting  the  complex  nature  of  the  volume  growth  which  depends  on 
such  factors  as  cell  loss,  cell  supporting  matrix,  cell  death,  slow  cell  division  due  to 
shortage  of  nutrients  or  oxygen,  etc.  (file:  growth_bal  1 12_031202.xls).  Data  on  the 
tumor  volume  growth  was  accumulated  continuously  and  to  followed  up  for  longer  time. 
Figure  2B  plots  the  measured  tumor  volume  as  a  function  of  time  with  the  most  up-to- 
date  data.  The  volume  doubling  time  for  the  late  growth  phase  is  determined  to  be  3  days, 
(file:  XC.08  I-125.xls). 

Fitting  of  the  in  vivo  tumor  growth  measured  in  7/8/91  on  the  p  line  (p443,  p445, 
Fig.  3  and  4  respectively)  showed  a  slightly  different  tumor  volume  doubling  time.  The 
volume  doubling  time  in  the  late  grow  phase  was  2.5  and  2.9  days  for  the  P445  and  P443 
lines,  respectively.  Only  the  P445  line  exhibited  two  different  growth  phases.  The 
doubling  time  for  the  initial  and  late  phases  was  1.6  days  and  2.9  days,  respectively.  It 
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appears  that  there  was  not  enough  data  taken  for  the  tumor  growth  in  the  time  range  of 
initial  growth.  For  the  late  growth  phase,  the  volume  doubling  time  obtained  from  this  set 
of  data  is  consistent  with  the  current  data,  (file:  growth_p_line_070891.xls). 

Therefore,  the  value  for  the  late  tumor  grow  volume-doubling  time  will  be 
determined  as  3.0  days.  The  initial  stage  volume  doubling  time  is  about  1.1  to  1.6  days, 
which  is  slightly  longer  than  the  in  vitro  cell  doubling  time  (20  hours).  Since  onset  of  in 
vivo  CLDRI  taken  place  at  about  19-21  days,  the  late  phase  volume-doubling  time  will 
be  used  for  the  model  calculation.  This  is  valid  if  the  tumor  density  is  assumed  to  be 
constant.  Cell  loss  and  death  and  reduced  cell  proliferation  due  to  short  supply  of 
nutrients  and  oxygen  will  slow  the  in  vitro  cell-doubling  rate.  The  growth  of  blood  vessel 
and  tumor  cell  supporting  matrix  complicates  the  relationship  between  the  cell  number 
and  tumor  size.  (Presumably,  the  outer  layer  may  grow  with  cells  of  in  vitro  doubling 
time,  but  the  inner  layer  which  has  less  sufficient  nutrient  may  have  cell  number  increase 
less  than  that  of  the  in  vitro  doubling  time).  Since  these  factors  are  not  taken  into  account 
explicitly  in  the  model  calculation,  the  use  of  the  late  phase  volume  doubling  time  as  the 
cell  number  doubling  time  the  model  calculation  is  a  reasonable  approximation,  as  further 
supported  by  the  data  described  next. 

c.  Tumor  cell  growth  based  on  in  vivo  tumor  growth  and  in  vitro  cell  assay 

Methods:  Tumor  cells  of  4500  to  7500  are  inoculated  to  the  head  of  the  rat.  Tumors  were 
allowed  to  grow  on  the  rat  until  19  to  21  days.  The  cells  from  the  tumor  at  19  to  21  days 
were  measured  by  cell  assay  of  the  tumor. 

Results:  1)  number  of  cells  measured  from  the  tumors  grown  on  the  rat’s  head  at  19-21 
days  spread  from  l-7xl06;  the  measured  cell  number  correlates  linearly  with  the 
measured  volume  of  the  tumor  (Fig.5).  Therefore,  the  use  of  late  phase  volume  doubling 
time  for  model  calculation  is  a  reasonable  one.  Because  the  volume  growth  as  described 
in  l.b,  progress  from  a  fast  early  to  a  late  slow  phase  within  this  period  of  time, 
assumption  of  a  simple  exponential  growth  from  day  of  tumor  cell  inoculation  to  the 
elapsed  days  of  19  -  21  day  will  not  be  valid.  If  this  assumption  was  used,  the  cell 
dumber  doubling  time  estimated  from  this  assay  technique  ranged  from  1.8  to  2.6  days. 
As  expected,  it  is  smaller  than  late  phase  doubling  time  of  3.0  days,  (file: 
cellnumber_at_rx  .xls) . 

2.  Sublethal  damage  repair  time  constant 

Methods:  Split  dose  experiment  was  performed  with  250  kV  x-rays  for  Ba-1 1 12  cell 
line.  A  total  dose  of  10  Gy  was  given.  The  survival  fraction  was  measured  as  a  function 
of  time  interval  between  two  5  Gy  irradiation.  The  cells  were  kept  at  22°C  during  the 
waiting  interval  to  prevent  proliferation.  According  to  incomplete  repair  model  we  have, 

ln[5(Q/5(0])  =1 
ln[S(°o)/S(0)]  6 
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where  S(t)  is  the  surviving  fraction  with  split  time  interval  t  and  |n  is  the  sub-lethal 
damage  repair  time  constant. 

Results:  By  fitting  the  measured  data  to  the  above  equation  (Fig.6),  the  in  vitro  repair 
half-time  for  the  BA1 1 12  cells  was  determined  to  be  20  minutes  (split_dose_bal  1 12.xls). 

3.  a  and  (3  for  the  linear-quadratic  cell  survival  model 

Methods:  The  parameters  of  a  and  (3  for  the  Ba-1 1 12  cell  line  for  1251,  103Pd,  and  250  kV 
photons  were  determined  from  the  survival  curves  obtained  under  acute  exposure 
condition  using  simulated  125I  and  103Pd  x-ray  beam  on  an  orthovoltage  unit. 

The  animals  chosen  for  the  in  vitro  irradiations  of  the  BA1 1 12  sarcoma  cells  were 
euthanasia  by  anesthetic  overdose  and  the  tumor  cells  were  then  removed  using  aseptic 
techniques.  A  single-cell  suspension  of  tumor  cells  was  suspended,  counted,  and  assayed 
for  viability  using  the  same  colony  formation  assay  used  for  cells  in  cultures.  The 
response  of  tumors  to  irradiation  will  be  assessed  by  the  ability  of  the  BA  1 1 12  tumor 
cells  to  form  colonies  in  cell  culture  after  in  vitro  irradiations  at  energies  equivalent  to 
125I  and  103Pd  under  aerobic  conditions,  and  leading  to  exponential  growth  at  different 
rates  (from  near  quiescence  to  full  exponential  growth  at  a  maximal  rate). 

To  obtain  the  a  and  (3,  the  acute  survival  curves  for  the  simulated  1251,  103Pd,  and 
the  250  kV  x-rays  were  fitted  to  the  LQ  cell  survival  equation.  Several  factors  need  t  be 
considered  for  such  a  fit.  1)  If  the  actual  cell  survival  fraction  is  used  in  the  fit,  the  least 
square  fit  will  automatically  favor  the  data  points  with  large  survival  fractions,  i.e. 
survival  data  at  low  doses  dominates  the  fitting  process.  2)  To  increase  the  importance  of 
large  dose  data  points,  the  logarithms  of  the  survival  fractions  were  used  to  fit  the  LQ 
model.  Although  it  gives  an  improved  fit,  it  does  tend  to  emphasize  the  large  dose  (small 
survival  fraction)  data,  of  which  have  usually  large  uncertainties.  3)  If  the  survival  curves 
of  1251,  103Pd,  and  the  250  kV  were  fitted  individually  to  LQ  equation,  the  a  and  (3  will 
assume  different  values  for  each  beam.  It  probably  represents  the  best  data  for  each  beam, 
and  the  use  of  these  parameters  in  model  calculations  for  each  beam  characteristics  may 
yield  the  best  agreements  with  experiments  for  the  given  beam.  However,  it  will  restrict 
the  model  calculation  applicable  for  only  the  beam  characteristics  that  has  the  fitted 
parameters.  The  model  will  loss  its  general  predictive  power. 

Initially,  we  have  taken  a  two-step  fitting  procedure.  In  the  first  step,  the 
logarithm  of  the  survival  data  was  fitted  to  the  LQ  equation  individually  for  each  beam. 
The  average  of  the  (3  from  these  fits  is  then  used  as  a  fixed  parameter  for  a  second  fit  of 
log  survival  curves  to  determine  the  a  for  each  beam.  The  rationale  was  that  beta 
represents  the  cell  death  due  to  sublethal  damages  that  is  intrinsic  to  each  cell  line. 
However,  this  procedure  would  lump  all  dose  rate  effects  to  the  parameter  a.  As  a  result, 
it  produced  quite  different  o/(3  ratio  for  each  beam.  Since  <x/f3  represents  the  dose  at 
which  the  cell  damage  from  single  hit  equals  to  that  from  double  hits,  regardless  of  dose 
rates  (?),  oe/(3  should  be  kept  constant  during  the  fit.  The  new  procedure  developed  for  the 
parameter  fitting  is  as  follows.  The  log  survival  fraction  from  all  three  beams  were  used 
simultaneously  in  a  least  square  fit  in  which  the  a  for  each  beam  and  cc/ (3  are  treated  as 
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the  global  fitting  parameters.  The  measured  survival  curves  and  the  calculated  curves 
using  the  fitted  parameters  were  given  in  Fig.7A.  The  a  and  P  determined  by  this 
procedure  are  given  in  the  following  table.  The  parameters  determined  from  a  global  fit 
in  which  the  P  was  assumed  to  be  the  same  for  all  beams  is  also  given  for  a  comparison 
(Fig.  7B).  It  is  interesting  to  note  that  the  values  of  the  a  and  P  determined  from  these 
two  fitting  procedures  are  not  that  different.  A  literature  search  on  the  radiobiological 
parameters  on  the  BA1 112  cell  line  did  not  yield  previously  determined  parameters  for 
the  LQ  cell  survival  model.  Further  investigation  on  the  process  of  determine  these 
parameters  from  the  measured  survival  curves  at  different  conditions  is  planned. 


Radiobiological  Parameters  Determined  from  Acute  Exposure  Experiments 


Parameters  for  Ba-1112 


Beam 

a 

(Gy1) 

P 

(Gy2) 

cx/p 

(Gy) 

Tumor 

Doubling  Time 
(day) 

Sub-lethal  Damage 
Repair  Half-time 
(hr) 

250  kV  old 

0.25 

0.041 

6.1 

3.0 

0.33 

125I  equivalent 

0.26 

0.043 

6.1 

3.0 

0.33 

103Pd  equivalent 

0.29 

0.048 

6.1 

3.0 

0.33 

250  kV  old 

0.23 

0.044 

5.1 

3.0 

0.33 

125I  New 

0.25 

0.044 

5.6 

3.0 

0,33 

103Pd  New 

0.32 

0.044 

7.3 

3.0 

0.33 

Figure  2 
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APPENDIX  V 

IN  VITRO  STAINING  PROTOCOL  FOR  FLOW  CYTOMETRY 

L  GROW  CHO  CELLS  AS  A  MONOLAYER  @  37°C  IN  BME  PLUS  15%  FBS,  VITAMINS,  AND  P/S. 
GROW  BA- 11 12  FROM  PRIMARY  EXPLANT  IN  D-MEM  PLUS  15%  FBS,  P/S,  Na-PYRUVATE, 
VITAMINS. 

2.  LET  CELLS  OBTAIN  LOGARITHMETIC  GROWTH.  THEN  LABEL  CELLS  WITH  IUdR  FOR  24  HRS.  IUdR 
CONCENTRATION  FOR  LABELLING  IS  10’5M  &  104M  IUdR. 

PE=  NO  ANTIBODY  TREATMENT,  CELLS  W/O  IUdR 

CONTROL^  SECONDARY  ANTIBODY  TREATMENT  ONLY,  CELLS  W/O  IUdR 
EXPERIMENT AL=  PRIMARY/SECONDARY  ANTIBODY  TREATMENT,  IUDR -LABELLED  CELLS 

3.  AFTER  24  HRS  INCORPORATION  TIME,  TRYPSINIZE  CELLS,  THEN  COUNT  CELLS  AND  DISTRIBUTE  3x10 
CELLS/TUBE  INTO  5  MICROCENTRIFUGE  TUBES.  (3-EXP,  1-CON,  1-PE)  FIX  THE  CONTROL  AND  EXPERIMENTAL 
TUBES  WITH  0.5  ML  OF  70%  ETHANOL  FOR  30  MINUTES  ON  ICE  (4°C).  LEAVE  PE  TUBE  IN  Pi/NaCL. 

4.  REMOVE  CELLS  FROM  ETHANOL  (3xlQ6  CELLS).  RESUSPEND  INTO  1.0  ML  OF  1.5M  HCL  @  20°C  FOR  20  MINUTES. 

5.  WASH  CELLS  TWICE  WITH  0.5  ML  Pi/NaCL 

PRIMARY  ANTIBODY 

6.  RESUSPEND  IN  A  Pi/NaCL  SOLUTION  CONTAINING  .5%  TWEEN-20,  .5%  BOVINE  SERUM  ALBUMIN  (MONO  I),  AND 
A  1:50  DILUTION  OF  A  MONOCLONAL  ANTIBODY  (RABBIT  ANTI  5  IODOURIDINE)  (0.4  ML)  FOR  1HR  AT  ROOM 
TEMPERATURE  (20°C). 

REMEMBER:  DO  NOT  ADD  THIS  ANTIBODY  TO  THE  CONTROL  TUBES! 

7.  WASH  CELLS  TWICE  AGAIN  WITH  0.5  ML  Pi/NaCL. 

SECONDARY  ANTIBODY 

8.  RESUSPEND  INTO  A  Pi/NaCL  SOLUTION  CONTAINING  .5%  TWEEN-20,  1%  NEUTRAL  GOAT  SERUM  (MONO  II), 
AND  A  1:100  DILUTION  OF  FLUORESCEIN-LABELLED  GOAT  ANTI-RABBIT  IgG  (0.4  ML)  FOR  30  MINUTES  OF 
ROOM  TEMPERATURE  (20°C). 

REMEMBER.  ADD  THIS  ANTIBODY  TO  CONTROL  AND  EXPERIMENTAL  TUBES. 


9.  WASH  CELLS  AGAIN  WITH  0.5  ML  Pi/NaCL 


10.  CELLS  MAY  BE  USED  FOR  DETERMINATION  OF  IUDR  CONTENT  BY 

PROPIDIUM  IODIDE  STAINING  AND  MEASUREMENT  BY  FLOW  CYTOMETRY  @ 

600nM  (RED  FLUORESENCE),  OR  LUMINESCENSCES  MEASUREMENTS  @  5 14nM(GREEN  FLUORESCENCE). 


In  Vivo  BrdU  Labeling  Protocol  For  Flow  Cytometric  Analysis 

Fixation  of  Cells  in  75%  Alcohol  (Step  1): 

An  injection  of  50  mg  per  kg  BrdU  is  given  Intraperitoneally,  and  at  appropriate  times  thereafter  a  rat  was  sacrificed  and  the  tumor  was 
excised,  cut  into  small  pieces,  washed  with  Hanks  basic  salt  solution,  and  incubated  at  37®C  for  10  minutes  in  Hanks  solution  containing 
0.05%  trypsin. 


The  resulting  cell  suspension  was  filtered  through  a  22  nylon  mesh  filter  and  centrifuge  at  1000  rpm  for  5  minutes. 
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Remove  the  supernatant  without  disturbing  the  pellet. 

FOLLOWING  STEPS  DONE  ON  ICE  (4°C) 

Add  1 .0ml  of  PBS/EDTA  solution  and  gently  stir  with  the  pipetman  in  order  to  get  the  single  cell  suspension. 

Add  4,0ml  of  PBS/EDTA 
Centrifuge  at  lOOOrpm  for  5  minutes. 

Remove  supernatant  again  without  disturbing  the  pellet 

Add  1.0ml  of  PBS/EDTA  solution  and  gently  stir  with  the  pipetman  in  order  to  get  the  single  cell  suspension. 

Now  add  1.0ml  of  95%  ethanol  drop  by  drop,  very  slowly  while  vortexing  the  pellet  at  a  very  low  speed,  (knob  setting  =  1) 

Now  add  1.0ml  of  95%  ethanol  drop  by  drop,  very  slowly  while  vortexing  the  pellet  at  a  very  low  speed,  (knob  setting  =  2) 

Now  add  1,0ml  of  95%  ethanol  drop  by  drop,  very  slowly  while  vortexing  the  pellet  at  a  very  low  speed,  (knob  setting  =  3) 

SUMMARY  NOTE:  3.0ml  of  95%  ethanol  should  be  added  in  three  steps,  1,0ml  at  each  step  drop  by  drop  while  vortexing  the  pellet. 
The  only  difference  is  the  increase  in  the  vortex  speed  from  1  to  3. 

The  tube  has  now  4,0ml  of  volume  (1 .0ml  of  PBS/EDTA  and  3,0ml  of  95%  ethanol,  (final  concentration  of  alcohol  is  75%) 
REMEMBER:  All  the  steps  above  must  be  done  on  ice,  including  putting  the  alcohol  on  ice  before  starting. 

NOTE:  You  can  keep  the  cells  at  4°C  for  one  month  until  proceed  to  the  next  step. 

Neutralization  of  Cells  (Step  2): 

Centrifuge  the  cells  @  1500  rpm  for  10  minutes. 

Remove  the  supernatant  without  disturbing  the  pellet. 

Add  1,0ml  PBS/EDTA  slowly  drop  by  drop  while  vortexing  the  pellet  gently,  (Do  not  use  the  pipetman  here  because  you  will  lose 
cells.) 

Add  4.0  ml  of  PBS/EDTA  to  bring  final  volume  to  5,0  ml. 

Centrifuge  at  1500  rpm  for  10  minutes. 

Remove  supernatant  with  a  Pasteur  pipette. 

Vortex  the  pellet. 

Add  1.0ml  of  2N  HCL /  0,5%TritonX-100  drop  by  drop  while  vortexing  at  low  speeds. 

Let  sample  sit  at  room  temperature  (20-25°C)  for  30  minutes. 

Centrifuge  at  1500  rpm  for  10  minutes. 

Remove  supernatant. 

Add  1.0ml  of  0,1  M  sodium  tetraborate  (Na2B407,10H20)  ph  8,5  drop  by  drop  while  vortexing  pellet  at  setting  2, 

Centrifuge  @  1500  rpm  for  10  minutes. 

Remove  supernatant. 

Add  1.0ml  of  PBS/0,5  %Tween-20/l%BS  A  solution  drop  by  drop  whiled  vortexing  the  cell  pellet, 

NOTE:  You  can  keep  the  cells  at  4°C  until  the  next  day. 
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Staining  The  Cells  with  Anti-BrdU  Antibody  (Step  31: 

Centrifuge  @  1500  rpm  for  10  minutes 
Remove  supernatant. 

Add  200  ul  of  PBS/0,5%Tween-20/l  %BSA  solution.  The  tumor  cells  are  rinsed  three  times  with  this  solution  before  incubation  with 
primary  antibody. 

Centrifuge  @  1500  rpm  for  10  minutes  after  last  rinse 

Add  0.5ml  of  the  primary  antibody  to  BrdU  (rat  monoclonal  Anti-BrdU  antibody  1:100  dilution  with  PBS/,05%Tween-20/l%  BSA 
solution). 

Incubate  for  60-90  minutes  at  room  temperature,  25°C,  with  a  gentle  tapping  of  the  tubes  every  15  minutes  to  get  uniform  staining  of  all 
cells.) 

Wash  cells  twice  again  with  0.5  ml  of  PBS/0.5%Tween-20/l%BSA  solution. 

Centrifuge  @  1500  rpm  for  10  minutes. 

Remove  supernatant. 

For  the  secondary  antibody  incubation,  resuspend  the  0.5ml  of  the  fluorochrome-conjugated  antibody,  (goat  anti-rat  antibody,  1 :200 
dilution  in  PBS/.05%Tween-20/l%  BSA  solution)  at  room  temperature,  25°C  for  45  minutes. 

Wash  cells  again  with  0.5  ml  of  PBS/.05%Tween-20/l%  BSA  solution. 

Propidium  Iodide  Staining-  Cellular  DNA  Content  (Step  41; 

Centrifuge,  decant  or  aspirate  fixative.  Insure  all  ethanol  has  been  removed  leaving  tube  upside  down  o  a  paper  towel 
Resuspend  cells  in  RNAse  at  2xl06  cells  or  minimum  of  0.4  ml.  Incubate  at  37°C  for  30  miutes. 

Add  Propidium  Iodide  for  a  final  concentration  of  0.05  mg/ml.  Keep  at  4°C  for  a  least  45  minutes  prior  to  flow  cytometric  analysis. 
Filter  samples  thorough  a  35  micron  nylon  mesh  immediately  prior  to  analysis. 
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Abstract 


An  experimental  technique  for  radiobiology  studies  of  brachytherapy  continuous  low 
dose  rate  irradiation  (CLDRI)  using  a  rat  solid  tumor  model  has  been  developed.  BA1 112 
tumors  were  grown  between  the  ears  of  the  14-week  old  male  WAG/rij  rats  by 
interdermal  inoculation.  Irradiation  of  the  tumor,  in  vivo,  by  the  125I  and  103Pd 
brachytherapy  sources  was  accomplished  by  using  source  afterloading  system,  which 
consists  of  a  lightweight  helmet  sutured  to  the  rat  and  a  nine-source  afterloading 
applicator.  The  applicator,  made  of  clear  polystyrene  with  loading  ports  for  nine  sources, 
has  a  12  mm  x  12  mm  opening  in  the  center  to  accommodate  the  tumor  and  its  growth 
during  the  irradiation  (the  diameter  of  a  typical  BA1 112  tumor  is  about  6  mm  when 
radiation  is  applied).  The  spatial  location  of  the  nine  sources  were  optimized  to  produce 
an  as  uniform  as  possible  three-dimensional  dose  distribution  to  the  central  portion  of  the 
applicator  for  both  12T  and  103Pd  sources.  Absolute  dose  delivered  by  the  applicator  was 
verified  by  point  dose  measurements  using  calibrated  TLD  in  a  polystyrene  phantom  that 
mimics  the  scattering  environment  of  the  tumor  on  the  rat.  This  experimental  technique 
was  applied  to  study  the  in  vivo  response  of  the  BA1 1 12  tumors  to  125I  and  103Pd 
treatment,  by  measuring  the  tumor  size  twice  weekly  during  the  treatment  until  each 
tumor  has  reached  a  maximum  volume  of  1  cm3  or  until  the  tumor  has  regressed  and  the 
animal  has  been  free  of  tumor  for  100  days.  The  technique  was  also  applied  to  for  tumor 
cure  studies  under  CLDRI  by  measuring  the  cell  survival  in  vitro  following  in  vivo 
irradiation  of  the  tumor.  In  this  application,  BA1 1 12  tumors  were  irradiated  at  a  low  dose 
rate  in  vivo  to  graded  doses  from  2  to  20  Gy,  the  tumours  were  then  removed  and  their 
colony  formation  ability  was  measured  using  an  in  vitro  assay  techniques,  (explicit 
results)  Description  of  the  technique  and  its  application  for  brachytherapy  radiobiological 
studies  are  presented. 
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length)  x  2.15  cm  (side-side  width)  x  2.2  cm  (height).  Priori  to  irradiation,  the  helmet  was 
sutured  to  the  rat’s  head  by  four  stitches  through  the  cartilage  of  the  ears  and  two  more 
stitches  just  behind  the  head  at  the  neck  (Fig.l).  A  seventh  suture  placed  under  the  tumor 
will  be  tied  to  the  central  bar  that  across  the  top  of  the  helmet,  thus  ensuring  the  tumor  is 
pulled  up  into  the  center  of  the  treatment  volume.  Radioactive  sources  of  125I  or  103Pd 
were  loaded  into  the  afterloading  applicator  in  the  hot  lab  and  the  applicator  was  then 
transported  to  the  animal  facility  to  be  afterloaded  into  the  helmet. 

Because  the  CLDRI  experiment  lasts  over  a  long  period  of  time,  special  radiation 
shielded  animal  housing  was  built  with  transparent  Pb-plastic  door,  which  provides 
adequate  light  and  air  for  the  rats  during  the  long  irradiation.  The  experimental  protocol 
was  reviewed  and  approved  by  the  Yale  University  Animal  Care  and  Use  Committee 
priori  to  the  start  of  the  experiments.  The  attachment  of  the  helmets  and  gesturing  during 
the  long  irradiations  has  shown  no  ill  effects  on  the  animals.  The  rats  have  gained  weight 
as  expected  and  have  gone  on  their  daily  activities. 

BA  1112  tumor  model  characteristics 

Tumors  used  for  in  vivo  tumor  cure  studies  have  been  implanted  by  the 
subcutaneous  inoculation  of  about  5000  tumor  cells  into  the  subcutaneous  tissues 
between  the  ears  of  14-week-old  male  WAG/rij  rats.  The  cells  were  obtained  from  a 
single  cell  suspension  of  BA1 1 12  cells  harvested  from  the  BA1 1 12  tumor  growing  on  the 
head  of  a  previously  inoculated  rat  at  an  elapsed  time  of  about  21  days.  The  rats  are  bred 
and  maintained  in  our  breeding  colony  under  SPF  conditions.  (Brief  description  (or 
reference)  of  the  BA1112  tumor  cell  line  and  its  characteristics  and  suitability  for 
brachy therapy  radiobiology  studies). 

Figure  3  plots  the  measured  tumor  volume  as  a  function  of  time.  The  tumor 
volume  growth  exhibited  two  growth  phases:  an  initial  fast  growing  phase  and  a  late  slow 
phase.  From  an  elapsed  time  of  15  days  on,  the  tumor  volume  grow  exponentially  with  an 
apparent  volume  doubling  time  of  about  3  days.  The  late  volume  growth  rate  is 
significantly  slower  than  the  in  vitro  tumor  cell  doubling  time  (about  20  hours),  reflecting 
the  complex  nature  of  the  volume  growth  that  is  affected  by  such  factors  as  cell  loss,  cell 
death,  prolonged  cell  cycle  due  to  local  environment  as  well  as  the  supply  of  nutrients 
from  the  tumor  cell  supporting  matrix. 

Tumor  cure  studies  under  125I  and  103Pd  CLDRI 

Two  techniques  have  been  developed  to  study  the  tumor  response  to 
brachytherapy  CLDRI  treatment  using  12SI  and  103Pd  sources.  In  the  first  approach, 
irradiation  of  the  tumor  was  applied  at  about  day  21  since  the  inoculation  of  the  tumor 
cells  to  the  rat.  Typical  tumor  volume  at  onset  of  radiation  is  about  100-200  mm3.  The 
response  of  the  tumors  to  CLDRI  was  analyzed  by  measuring  the  tumors  twice  weekly 
during  the  irradiation  until  each  tumor  has  reached  a  maximum  volume  of  1  cm3  or  until 
the  tumor  has  regressed  and  the  animal  has  been  free  of  tumor  for  100  days.  The  pattern 
of  tumor  growth  was  analyzed  and  TCD50  were  calculated  for  comparing  the  effects  of 
different  treatments. 
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In  the  second  approach,  an  in  vitro  tumor  cell  survival  assay  from  tumor  irradiated 
in  vivo  was  used.  Irradiation  of  the  tumors  growing  on  the  rats,  to  different  radiation 
dose,  commenced  at  approximately  3  weeks  after  inoculation.  The  rats  were  killed  at  the 
end  of  the  irradiation  and  the  tumor  cells  were  suspended,  counted,  and  assayed  for 
viability  using  the  same  colony  formation  assay  used  for  cells  from  cultures.  Analyses  of 
cell  yield  were  performed  to  allow  for  the  loss  of  cells  during  the  protracted  irradiations. 
This  endpoint  was  used  to  measure  the  cell  survival  curves  after  relatively  short,  graded 
treatment  times  (hours  to  days). 

Experiments  were  conducted  using  125I  sources  with  mean  initial  dose  rate  of  8 
cGy/hr  and  16  cGy/hr,  respectively.  For  1®Pd  sources,  a  mean  initial  dose  rate  of  16 
cGy/hr  was  used. 


Results  and  Discussion 

The  intrinsic  initial  iso-dose-rate  distributions  (in  cGy/hr)  in  the  planes  parallel  to 
the  base  of  the  applicator  for  sources  of  unit  air-kerma  strength  are  shown  in  Fig.  4  for 
125I  seeds  (Draximage  model  LS-1)  and  in  Fig.  5  for  103Pd  (Theragenic  Model  200)  seeds. 
Eight  representative  planes  parallel  to  the  base  of  the  applicator,  located  at  -1mm  (A), 
0mm  (B),  1mm  (C),  2mm  (D),  3mm  (E),  4mm  (F),  5mm  (G),  and  6mm  (H)  above  the 
base,  are  shown  in  each  graph.  Figure  6  plots  the  cumulative  dose-volume  histogram  for  a 
cylindrical  volume  at  the  center  of  the  applicator  (the  cylinder  has  diameter  of  10  mm  and 
a  height  of  10  mm).  It  is  seen  that  the  iso-dose-rate  distribution  inside  the  open  portion  of 
the  application  has  a  uniformity  between  95%  and  150%  for  both  125I  and  1<33Pd  seeds. 
Table  I  compares  the  dose  measured  in  polystyrene  to  the  dose  calculated  to  water  at  the 
same  measurement  point  for  six  applicators.  The  measured  dose  in  polystyrene  is  on  the 
order  of  5%  higher  than  the  calculated  dose. 

Figure  7  and  8  show  tumor  growth  under  CLDRI  for  125I  with  initial  dose  rate  of  8 
cGy/hr  and  16  cGy/hr,  respectively.  At  the  initial  dose  rate  of  8  cGy/hr,  the  tumors  were 
not  cured.  However,  the  tumor  growth  was  slowed  from  a  volume-doubling  time  of  3 
days  to  an  apparent  volume-doubling  time  of  21  days  under  CLDRI.  When  the  initial 
dose  rate  is  increased  to  16  cGy/hr,  tumor  cure  was  observed  for  the  initial  four  animals. 
The  experiments  on  four  more  animals  are  being  conducted  which  will  complete  by  the 
end  of  June,  2003.  The  prelimary  data  indicates  that  there  is  a  75%  tumor  cure  for  the 
BA1112  tumor  when  irradiated  by  125I  with  initial  dose  rate  of  16  cGy. 

Figure  9  plots  the  similar  in  vivo  tumor  cure  study  for  the  103Pd  sources  with 
initial  dose  rate  of  16  cGy.  The  tumor  cure  rate  is  50%  based  on  the  initial  experiments. 

Figure  10  plots  the  tumor  growth  as  a  function  of  treatment  time  when  the  tumor 
was  irradiated  by  a  mixture  of  125I  and  103Pd  sources.  The  total  initial  dose  rate  was  16 
cGy/hr  with  125I  and  103Pd  each  contributing  8  cGy/hr.  Under  this  condition  of  CLDRI, 
the  tumor  cure  rate  was  75%  from  the  prelimary  data. 
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Figure  3 
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Figure  6 
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